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I. INTRODUCTION 

Although the Standard Model (SM) has successfully accounted for all low energy CP- 
violating phenomena thus far observed, it is insufficient as the source of CP- violation needed 
to explain the cosmological baryon asymmetry in our universe According to Sakharov 
one of the conditions necessary to generate such an asymmetry is that the elementary inter- 
action violates charge conjugation symmetry (C) and the combined CP symmetry (where 
P is the parity symmetry). However, the size of the asymmetry needed for this cannot be 
derived in model calculations based on the SM and new sources of CP violation have 
been sought for many years in particle physics experiments. Prominent among these are the 
rare decays K — > -nvv which are sensitive to new physics involving both CP-violating and 
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CP-conserving interactions. In this paper, we present a detailed description of the previ- 
ously reported measurement of the reaction K + — > tt + vv performed by Brookhaven National 
Laboratory (BNL) experiment BNL-E949 4J. This paper is arranged as follows. We first 
briefly review CP violation and rare kaon decays, with an emphasis on K + — > -k + vv de- 
cays. We then describe previous results on this reaction and discuss the sources of potential 
background and the methods for suppressing backgrounds. We also discuss the design of 
the K + beamline, the detector and the selection criteria used in data analysis and describe 
the methods used for estimating background levels and for evaluating the acceptance. After 
examining the signal region, we present the method used for extracting the branching ratio, 
making full use of our knowledge of the background in the signal region. In the last section, 
we show how the measurement of B(K + — > tt + uu) impacts the search for new physics beyond 
the SM. 



A. CP Violation and the Rare Decay K + — ► tt + vv 



Standard Model CP violation arises from a complex phase in the three-generation quark 
mixing matrix 5]. In the Wolfenstein parameterization |f| of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, the parameters can be written in powers of A = sin^c ~ 0.22: 

( T/ - T/ T/ ^ ( l-X 2 /2 X A\ 3 {p-i v )^ 

-A l-A 2 /2 AA 2 (1) 

A\ 3 {l-p-iT]) -AX 2 1 

Where A, A, p and 77 are real numbers. CP invariance of the Lagrangian for weak interactions 
is violated when the CKM matrix is complex. The parameter 77 quantifies CP violation in 
the SM. 

The unitarity of the CKM matrix implies six unitarity conditions, which can be repre- 
sented graphically in the form of triangles, all of which must have the same area. The area of 
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these triangles is equal to one half of the Jarlskog invariant, Jqp 
property VW = 1 to the CKM matrix in ([T]) yields 



Applying the unitarity 



v: b v ud + v; h v cd + v; h v td ~ vi - xvx + v t 



<cb v cd 



tb v td 



ub 



cb 



'td 



0, (2) 



where the approximations V ud ~ V£ ~ 1 and V cd ~ —A have been made. This equation 
can be represented graphically, as shown in Fig. [fl where we have divided all of the sides 



(P.T1) 




(0,0) 



(1,0) 



FIG. 1: Unitarity triangles in the p — fj plane. Two sides of the triangle can be expressed by the 
CKM matrix elements |I4d|/^4A 3 and \V u b/V c i\/\, respectively, where A and A are parameters in 
the Wolfenstein parameterization. 

by W* b . The apex of the triangle is given by two Wolfenstein parameters, p and fj, where 



B's and K's are so far the only two mesons showing evidence of CP violation in their 
decay processes. Whether or not the observed CP violation can be completely explained by 
the CKM phase within the SM can be probed by the independent determination of p and rj, 
from B and K decays as shown in Fig. [2j Two sensitive methods for making the comparison 
are: 

• A comparison of angle /3 from the ratio B(K L — > tt°uu)/B(K + — > ir + uu) with that 
from the CP violating asymmetry («4.Qp) in the decay B^ — > J/ipK®; and, 

• A comparison of the magnitude \Vtd\ from K + — > tt + vv with that from the mixing 
frequencies of B s and B^ mesons, expressed in terms of the ratio of the mass differences, 



Although the decay K + — > tt + uu is a flavor changing neutral current (FCNC) process 
prohibited at tree level in the SM, it is allowed at the one-loop level. In leading order, it 
is described by a "Box" diagram and two "Z-penguin" diagrams, as shown in Fig. [31 The 
weak amplitude for this process is represented as 



p = p(l- A 2 /2) and fj = rj(l - A 2 



/2) i. 



AM B JAM Bi 




(3) 



S 



(p/n) 



K + -Mt + VV 




(0,0) 



(1,0) (p,„0) 



FIG. 2: Unitarity triangle determined by B and K decays. The parameters p and f) can be 
determined in two ways: the angle /3 from the CP violating asymmetry in the decay B® — > J/ipK®, 
and from the length of the side from AMbJAM^ in B — B mixing; the height of the triangle 
from B(K^ — > tx^vv) and the radius of a circle centered at (po>0) from B(K + — > 7r + z/P). 



where V^-'s are the CKM matrix elements, 7 M 's are the Dirac matrices, is the momentum 
transfer, and m^'s are quark masses. M. vanishes if all of the quark masses, mj, are equal, 
because of the unitarity of the CKM matrix. However, the breaking of flavor symmetry, 
which results in the variation of quark masses, allows this decay to proceed at a very small 
rate. The top quark provides the dominant contribution to the K + — > 7r + uu branching ratio 
due to its very large mass in spite of the small coupling of top to down quarks (Vtd) in the 
CKM matrix. 

Following Ref. [sj], the branching ratio for K + — > -k + vv is calculated as follows. The 
effective Hamiltonian can be written in the SM as 



n/SM 
H eff 



a 



a/2 2tt sin 2 9^ , 

l=e,fj,,T 

in next-to-leading order (NLO), where 



KVciXL + v*v td x(x t ))(sd) v - A (m)v~A, 



(4) 



and 



X( Xt ) EE X ( Xt ) + ^Ijflt) « Tlx ■ X ( Xt ) 



X (x t ) = Co(x t ) - AB {x t ), Vx = 0.995. 



(5) 



(6) 



Bo(xj) and Co(xj) in ([6]) are functions of Xj = m 2 /M^, and were derived for the first time 
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u,c,t 




V ^- V 

FIG. 3: Second-order weak processes that contribute to the K + — ► ix^vv branching ratio: the 
"Box" diagram (upper) and two "Z-penguin" diagrams (bottom). 



by Inami and Lim in 1981 [10]. The coefficient X^ L and the function X(xA are the charm 
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and top quark contributions, including QCD corrections at NNLO 

With the top quark mass in the minimal subtraction scheme m t (m t 
2.2) GeV 3, 

X(x t ) = 1.464 ±0.025 
is obtained. 

The perturbative charm contribution gives the largest theoretical uncertainty and can be 
described in terms of the parameter 



'162.3 ± 



(7) 



Pc(X) = ^X^ + ^ h ] 



0.34 ±0.04, 



(8) 



where the error is obtained by varying the charm mass, m c , the scale factor, \x c = 0(m c 
and the coupling constant, a s (M|), by reasonable amounts. One obtains 



where 



7T vv ) = Kj 



ImA* 



V (P c (X) + SP c , u ) + ^X( Xt ] 



A 



3a 2 B(K + -> 7T e+u) . , „ 
k+ = r+ _ v „ . t - ^A 8 = (5.26 ± 0.08) x 10" 11 



2vr 2 sin 4 6 



w 



A 



0.2257 



(9) 



(10) 
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SP c ,u = 0.04 ±0.02 comprises the long-distance contribution calculated in Ref. [if]], and the 
Aj's (= VjgVjd) are from the CKM matrix elements. The r + (= 0.901) represents isospin 



breaking corrections in relating K + — > n + vv to the well-measured decay K + — > 7r°e + i> 17]. 
In obtaining the numerical value in fTTUl) . we used [3] 

sin 2 Q w = 0.231, a = ——, B(K + n°e + u e ) = (4.98 ± 0.07) x 10~ 2 . (11) 

127.9 

Expression (Q describes in the p — fj plane an ellipse with a small eccentricity, namely 

W + i»-M'= ikW ' (12) 

where 

"°= 1+ | W X H ' ^^-Tj • K+ =A?' (13) 
Using d^J) and varying m t , |V c b|, -P c (^) and |Vtd|, which is constrained by IV^/V^I and 
B — B mixing in the p — fj plane, the branching ratio of K + — > 7r + z/P is predicted to be 

B(K + -> vr+^z/) = (0.74 ± 0.20) x 10" 10 (14) 

within the SM. It should be noted that, of the uncertainty of 27% in (JHJ), the theoretical 
uncertainty is ~6% at present, mainly due to the uncertainty in the charm quark mass. 

Theoretically a precise measurement of B(K + — ► tt + vv) is one of the cleanest ways to 
extract \V t d\- This is due to the following factors: 



the long-distance contributions to the branching ratio are small 19] and under control 



the most recent calculation gives a contribution of (+6±3)% to the branching ratio 



3; 



the uncertainty from the hadronic matrix element has been reduced to <1% by recent 

□ 

theoretical and experimental developments [20(, and 



HE? 



• the recent NNLO calculation 

mm 

has reduced the total theoretical uncertainties to 
~6%, i.e. relatively small and reliably calculated as compared with the uncertainties 
present in other K and B decays. 

If a precise measurement of the neutral analog K\ — > -n^vv could also be made, the intrinsic 
theoretical error on \V t d\ could be reduced to ~1% [is! ]. 

As determinations of 5-system parameters become increasingly precise, the uncertainty 
on the SM prediction for K + — > 7r + z/P will approach the current theoretical accuracy of 

11 



~6%. A correspondingly precise measurement of the K + — * tx + vv branching ratio there- 
fore provides a stringent test of the SM and probes for new physics. There have been 
numerous predictions for K + — > tx + vv in the models beyond the SM and applications of 
the measured branching ratio to constrain new models. These include the Minimal Su- 
persymmetric Standard Model with 



221 ] and without [221, |23| new sources of flavor- or 
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CP-violation, generic SUSY with minimal particle content [24L SUSY with non-universal 



A terms 
color models 
toquarks 



29 



251 ] . SUSY with broken R-parity [26|, |27j, top color [281 ] . top color-assisted techni 



30|, multiscale walking technicolor [31 L f° ur generation models 32|, lep- 



33], Left-Right model with right-handed Z' [34j, extension of the SM to a gauge 

a multi-Higgs multiplet mod el 361 ] . light sgoldstinos [37J, 



theory with J = mesons 
universal extra dimensions 



nano 



5-dimensional split fermions 



39] . a Randell-Sundrum sce- 



40J, a littlest Higgs model 4l|, |42J, non-standard neutrino interactions 43], and a 



minimal 3-3-1 model 



44). 



B. History of K + — ► 7r + uu Experiments 

Searches for this process which began over 35 years ago have used stopped- K + beams. It 
was believed at the time of the first of these that the branching ratio might be as high as a 



few xl0~ 5 45]. It was recognized that even at this level, a poor-signature process such as 
K + — > 7t + uj? would need effective particle identification, precise kinematic measurement and 
the ability to veto extra charged and neutral tracks to discriminate it from common decay 
modes such as K + — > {i^v^ and K + — > 7r + 7r° (referred to as K^ 2 and K %2 , respectively). 



The earliest published result was from a heavy liquid bubble chamber experiment [46] at the 
Argonne Zero Gradient Synchrotron, in which a 90% CL upper limit B(K + — > -k^vv) < 10~ 4 
was obtained. In that paper it was recognized that K n2 decay in flight and hadronic 7r + 
interaction in the detector were dangerous sources of potential background. 

The final analysis of the Argonne experiment improved the limit to 5.7 x 10 -5 \^\ . but 
before it appeared in print, a subsequent counter/spark-chamber experiment at the Berkeley 



48] . However this experiment was sensitive to only 



Bevatron improved the limit to 1.4 x 10~ 6 
the most energetic of ir + , whereas the bubble chamber experiment covered a wide kinematic 
range. In addition to the background from common K + decay modes, this experiment 
considered possible background from K + charge exchange in the stopping target followed 
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by K\ — > ir + e~u e , and from beam ir + which scattered into the detector. The Chicago- 
Berkeley group continued their program with a setup sensitive to tt + in the kinetic energy 
range 60-105 MeV, i.e. below that of the potential background process K n2 rather than 
above it. This required reconfiguring their photon veto system so that it became nearly 
hermetic. Combining results from the two configurations, the branching ratio upper limit 



was improved slightly to 5.6 x 10 



49). 



About a decade later, an experiment at the KEK Proton Synchrotron improved the limit 
to 1.4 x 10 -7 j|3(3]. The technique of waveform digitization to record the ir + — > fi + — > e + 
decay chain was introduced for the first time. This experiment was sensitive only to the 7r + 
with momenta greater than that from K n2 (referred to as the "7ri/z/(l)" region) and its setup 



resembled that of Ref. 



48). 



The BNL series of experiments was initiated in the early 1980 's. They were based on a 
large-acceptance solenoidal spectrometer with a hermetic photon veto situated at the end 
of a highly pure, very intense stopped-.fr + beam 5l| from the BNL Alternating Gradient 
Synchrotron (AGS). The experimental signature of the K + — > 7i + uu decay was a single ir + 
track with tt + momentum less than 227 MeV/c plus no other particle from a K + decay. 
Fig. H] shows momentum spectra of major decay modes of K + . 



The first p 
of 2.4 x 10~ 9 
CL 



rase of E787 in 1988-91 achieved a 90% CL upper limit on the branching ratio 



521 ]. using data from the t[vu{1) region. A separate limit of 1.7 x 10 8 at 90% 



53] was extracted from the kinematic region in which the ir + is softer than that of the ir~ 



from K n 2 (referred to as the "7rz/i/(2)" region). This program completed the identification of 
backgrounds needed to reach the 10~ 10 level of sensitivity and developed methods to reliably 
measure them. 

A major upgrade of both the beam line and the detector was undertaken between 1992 
and 1994. The search for K + — > tt + uu resumed in 1995 and continued through 1998. The 
limit on the branching ratio from the 7rz^(2) region was improved by an order of magnitude 
to 2.2 x 10~ 9 at 90% CL 54j], but the major output of this series of runs was the observation of 
two clean K + — > tx + vv events 55| in the irvvil) region and a measurement of the branching 
ratio B(K+ -> ti + vv) = {1.57toll) x 10 ~ 10 - The BNL-E787 detector was upgraded again 
over the period from 1999-2001. The E949 experiment was proposed to use this detector to 
run for 60 weeks. After the first 12 weeks of running in 2002 no further funds were provided 
to complete the experiment. Based on the collected BNL-E949 data, the first result was 
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FIG. 4: Momentum spectra (in MeV/c) of charged particles from K + decays in the rest frame. 
The values in the parentheses represent the branching ratios of the decay modes . The hatched 
spectrum shows the ir + momentum from K + — > tt^vv decay assuming the V — A interaction. 

already published in 2004 [4]. This paper provides an extended and detailed description of 
the detector and data analysis techniques used to produce the E949 result. 



II. EXPERIMENTAL METHOD 



A. Overview 



E949 (BNL-E949) which succeeded BNL-E787 had a sensitivity goal of detecting ten SM 



signal events 56|. E949 employed a low momentum beam of K +, s which were degraded 
and stopped in the detector. Measurement of the K + — > tt + vv decay involved observation 
of the daughter ir + in the absence of other coincident activity. The tt + was identified by 
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its kinematic features obtained from energy, momentum and range measurements, and by 
the observation of a ir + — > fi + — > e + decay sequence. Since the signal was expected at the 
10~ 10 level, the detector was designed to have powerful tt + identification for rejection of 
K^ 2 and K + — > yU + ^ M 7 decays (K^), 4-7r solid angle photon detection coverage for vetoing 
K V 2 decays, and efficient K + identification system for eliminating beam-related backgrounds. 

The entire E949 spectrometer was surrounded by a 1 Tesla solenoidal magnetic field along 
the beam direction. The coordinate of detector used a Cartesian coordinate system in which 
the origin was at the center of the target; the +z axis was along the incident beam direction 
and the +y axis in the vertical up direction as shown in Fig. [5j Under this coordinate system, 
the azimuthal angle of a track was defined as the arctangent ofy/x and, the polar angle 9 



was defined as the an g! 
elsewhere |57 
2001 56 



58 



59 



e with the +z axis. Many detector components have been discussed 
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61 



62j. Fig. [5] shows the E949 detector after upgrades (1999- 



63j) with the following improved components: photon veto detection efficiency, 
tracking and trigger efficiency, and data acquisition (DAQ) live time. E949 was designed to 
run at the same instantaneous rate as E787, and to achieve a factor of five improvement in 
sensitivity, through the use of a higher duty factor and reduced K + momentum for a higher 
stopping fraction. The higher duty factor was not achieved in the engineering run in 2001 or 
the first physics run in 2002 due to a broken motor generator set that supplied power to the 
AGS. The regular supply was removed from operation on August 3, 2001 and the backup 



was used during the rest of 2001 and 2002. E949 ran at about twice the beam rate of E787. 



B. Accelerator and K + Beam Line 

The K + beam was produced by a high- intensity proton beam from the AGS at BNL: the 
entire AGS beam of 65xl0 12 protons (Tp/spill) at a momentum of 21.5 GeV/c was delivered 
to the E949 K + production target. Prior to 2001 the AGS typically ran at 24 GeV/c, but at 
this momentum the longest spill achievable was 0.5 sec. Combined with the longer cycle time 
(3.2 s between spills, as compared to 2.3 s), the duty factor at 24 GeV/c was unacceptably 
low. By lowering the proton momentum to 21.5 GeV/c, the spill length was increased to a 
maximum of 2.2 s, resulting in a duty factor of 2.2 s/5.4 s. At this lower proton momentum 
the production of 710 MeV/c K +, s was reduced by 10%. The K + production target was 
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(a) Barrel 
Veto 



Downstream 
Photon Veto 





BVL 

Barrel / Drift 

Veto / Chamber \ x l RSSC 
Range Stack Target 



FIG. 5: Schematic side (a) and end (b) views of the upper half of the E949 detector. Illustrated 
in this figure, an incoming K + that traverses all the beam instruments, stops in the target and 
undergoes the decay K + — > 7r + 7r°. The outgoing 7r + and one photon from tt° — ► 77 are also shown. 
The detector elements and acronyms are described in detail in the text. 

made of 2/3 of an interaction length of platinum (6 cm along the beam direction), and was 
located on a water-cooled copper base. At the typical AGS running condition 65 Tp on 
the production target per 2.2 s spill, the maximum target temperature was measured to be 
~700°C. 



The Low Energy Separated Beam 64j (LESB III) collected and transported K +, s emitted 
at 0° (along with 500 ir +, s and 500 protons per K + ), and momentum-selected by the first 
dipole magnet. Two electro- magneto-static separators swept 7r + 's and protons out of the K + 
beam axis. The resulting beam was further selected by a second dipole magnet. LESB III 
had a total length of 19.6 m from the production target to the E949 target with an angular 
acceptance of 12 msr and a momentum acceptance of 4.5% FWHM at a mean momentum of 
710 MeV/c. During most of the 2002 running period the first separator voltage was lowered 
from the standard voltage of 600 kV to ~250 kV due to high voltage discharges. Under 
these conditions a K + : ir + ratio in the beam of 3:1 was achieved with a 40% loss in K + flux 
(typically E787 ran with 4:1). Proton contamination was suppressed to a negligible level by 
the separators. At the same typical AGS running condition described previously, 1.3 x 10 7 
K + 's were transported through the beam line. 

The typical conditions during the 2002 run had 3.5xl0 6 K +, s entering the E949 target 
every spill. This corresponded to a rate of 1.6xl0 6 i^ + 's/sec during the 2.2 sec spill. The 
effective spill length was actually 8% shorter (2.0 sec) due to some residual modulation of 
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the beam intensity out of the AGS. The typical instantaneous rate of beam particles at the 
Cerenkov counter was 6.3 MHz of K + 's and 1.5-2.5 MHz of 7r + 's. 



C. Beam Counters 

The incoming K + beam traversed a scintillation counter (BO), a Cerenkov counter, two 
beam wire proportional chambers (BWPC's), a passive BeO degrader, an active degrader 
(AD) and a beam hodoscope (B4) as depicted in Fig. [5j The BWPC's and BO counter are 
not shown in Fig. [5J 

The BO counter, which was a 30.5-cm long, 0.6-cm thick and 7.6-cm wide Bicron BC408 
plastic scintillator, was located just downstream of the last quadrupole magnet and counted 
all charged particles in the beam. It was read out by an analog-to-digital converter (ADC), 
a time-to-digital converter (TDC) and a 500 MHz transient digitizers based on gallium- 



arsenide charge-coupled device (CCD) [57|. The Cerenkov counter [51j located just down- 
stream of the BO counter identified particles as K +1 s or 7r +, s. The Cerenkov light from the 
K + (ir + ) was transmitted (internally reflected) at the downstream surface of the Cerenkov 
radiator and read out with 14 U K Cerenkov" (C K ) and 14 "vr Cerenkov" (C n ) EMI9954KB 
PMT's. The PMT signals were split, with 90% sent to TDC's via fast LRS3412 discrimina- 
tors and 10% to a xlO amplifier. The amplifier output was sent to CCD's. The pulse-height 
information in every 2 ns interval was recorded to reproduce the time development of the 
pulses and to detect two particles close in time to each other. The multiplicity output of 
Cx(C n ) PMT's was discriminated (typical threshold: n > 5) to identify K + 's and 7r + 's in 



the trigger (KB and 7rB defined in Section Til Hj) . 



The two BWPC's were located downstream of the Cerenkov counter to monitor the 
beam profile and identify multiple incoming particles. The first chamber (BWPC1) was 
located 168.5 cm upstream of the target entrance and contained three planes of sense wires: 
vertical (x-plane) and ±45° to the vertical (u- and ^-planes). The sense wires were 12- 
/im-diameter gold-plated tungsten. The x-, u- and u-planes had 144, 60 and 60 readout 
channels, respectively, with a 1.27 mm wire spacing. In the u- and f-planes, pair of wires were 
multiplexed in one readout channel. The active area was 17.8 cm (horizontal) by 5.08 cm 
(vertical). The cathode foils were 25-/xm thick aluminized mylar coated with carbon. The 
anode-cathode distance was 3.18 mm, and the total thickness of BWPC1 was approximately 
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56 mm. The second chamber (BWPC2) was located 1.0 m downstream of BWPC1 and also 
contained three planes (x, u and v). The direction of the sense wires was vertical (rr-plane) 
and ±60° to the vertical (u- and -u-planes). Each plane had 120 active sense wires with a 
0.8- mm wire spacing. Among the 120 wires, the central 72 ones were multiplexed by 3 and 
the remaining were multiplexed by 6 in the readout channels, yielding a total of 32 readout 
channels for each plane. The cathode foils were 8-/xm single-sided aluminized mylar coated 
with carbon. The anode-cathode distance was 1.6 mm. Both chambers were filled with a 
recirculated mixture of CF4 (80%) and Isobutane (20%). 

Downstream of the BWPC's a degrader slowed the K +, s so that they stopped in the center 
of the scintillator fiber target. The upstream section of the degrader was inactive, consisting 
of 11.11 cm of beryllium oxide (BeO) and 4.76 mm of Lucite. The high density (3.0 g/cm 3 ) 
and low atomic number of BeO was used to minimize multiple scattering. The AD consisting 
of 40 layers of 2 mm thick disks of Bicron BC404 scintillator (13.9 cm diameter) alternating 
with 2.2-mm thick copper disks (13.6 cm diameter) was divided into 12 azimuthal segments 
with readout to a single Hamamatsu R1924 PMT through fourteen 1-mm-diameter Bicron 
BCF99-29-AA-MC wave length shifting (WLS) fibers. The PMT outputs were provided to 
TDC's, CCD's and a 4-fold analog sum that was provided to an ADC. These measurements 
enabled the AD to identify the beam particles and to detect activity coincident with K + 
decays. 

Downstream of the degrader the B4 hodoscope detected the entrance position of the 
incoming particle in the target and identified the particle type by measuring its energy loss. 
The B4 hodoscope consisted of two planes, u and v, with about a 11.8-cm diameter oriented 
at a ±33.50° angle with respect to the horizontal axis. Each plane had 16 Bicron BC404 
scintillator fingers with a 7.2-mm pitch. The cross section of each finger had a 'Z-shape' 
with a 6.4-mm thick middle part and 3.2-mm thick edges. This shape reduced inactive 
regions and improved the spatial resolution. Three Bicron BCF99-29-AA-MC WLS fibers 
were embedded in each finger and connected to a single Hamamatsu H3165-10 PMT that 
was read out by TDC's, ADC's and CCD's. At the same position as the B4 hodoscope, but 
at larger radius was an annular scintillator counter, the ring veto (RV). The RV was designed 
to veto particles that passed through perimeter of the B4 hodoscope. The RV was composed 
of two 180° arcs of 3.3 mm thick Bicron BC404 scintillator with an inner diameter varying 
from 11.9 cm to 12.0 cm and an outer diameter of 14.6 cm. The two RV elements were 
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readout by Hamamatsu H3165-10 PMT's and the signals were split three ways to ADC's, 
TDC's and CCD's. 



D. Target 

The target consisted of 413 Bicron BCF10 scintillating fibers of 5-mm square cross section 
and 3.1-m length that were bundled to form a 12-cm-diameter cylinder. A number of 1-mm, 
2-mm and 3. 5-mm square scintillating fibers (called "edge fibers") filled the gaps near the 
outer edge of the target. Each of the 5-mm fibers was connected to a Hamamatsu R1635-02 
PMT, whereas the adjacent edge fibers were grouped onto 16 PMT's, providing signal read 
out by ADC's, TDC's and CCD's. 

The fiducial region of the target was defined by two layers of six plastic-scintillating 
counters surrounding the target. The inner scintillators, called I-Counters (IC's), helped 
to define the fiducial volume and the logical OR of the six IC's (JC for trigger condition) 
used by the trigger for this purpose. The IC's were 6.4-mm thick (with an inner radius of 
6.0 cm) and extended 24 cm from the upstream face of the target. The outer scintillators, 
called V-Counters (VC's), overlapped the downstream edge of the IC's by 6 mm, and served 
to detect particles that decayed downstream of the fiducial region of the target. The VC's 
consisted of six 5-mm thick and 1.96-m long scintillators, and were staggered azimuthally 
with respect to the IC's. Each IC and VC element was instrumented with an EMI 9954KB 
PMT which was read out by an ADC, TDC and a 500 MHz transient digitizer (TD) based 
on a flash ADC 



Approximately 27% of the incident K +, s (typically 3.5 x 10 6 i^ + 's/spill) penetrated far 
enough into the target to satisfy the online target criteria for KB defined in Section III H 21 
The remaining K +, s either underwent decay-in-flight, nuclear interaction in the degrader 
or scattered in the material of the beam instrumentation and did not reach the target. It 
should be noted that some of the K + 's (< 25%) that satisfied the online KB requirement 
didn't stop in the target, and a factor for the stopping fraction of K +, s was introduced 
as described in Section [III H 71 The K +, s deposited an average energy of 100 MeV in the 
scintillating fiber target when coming to rest in the center of the target fiducial volume. The 
low velocity K +, s typically lost 5-40 MeV in each fiber, while the nearly minimum ionizing 
(MIP) 7r + 's from K + decays deposited about 1 MeV per fiber, as they passed transversely 
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through the fibers. 



E. Drift Chamber 



The drift chamber, called the "Ultra Thin Chamber" (UTC) [59j], was located just outside 
of the IC. The primary functions of the UTC were to measure the momenta of charged 
particles and to provide tracking between the target and the Range Stack (RS). The UTC 
had inner and outer radii of 7.85 cm and 43.31 cm, respectively. Twelve layers of 5-8 mm 
drift cells were grouped into three super-layers (each consisting of four layers) with active 
lengths of 38.8 cm (inner), 44.8 cm (middle) and 50.8 cm (outer). The super-layers were 
filled with a 49.8%:49.8%:0.4% mixture of argon, ethane and ethanol. Each anode wire was 
instrumented with an ADC and a TDC. The drift time to the anode wires was used to 
determine the (x, y) positions for the charged track. At the inner and outer radii of each 
super-layer were cathode foils, with 7 mm helical cathode strips at a ~ 45° pitch angle. Each 
cathode strip was instrumented with an ADC and a TDC. The combination of the charge 
centroid on the cathode strips and anode hits provided the z hit position. There were two 
inactive regions filled with nitrogen gas between the three super-layers. Differential pressures 
of ~2 mbar in the five gas volumes supported the cathode foils. The total mass of the UTC 
(excluding the inner and outer support tubes with the attached foils) amounted to 2xl0 -3 
radiation lengths. The UTC position resolutions were approximately 175 fim for x and y 
and 1 mm for z. 



F. Range Stack 

Located just outside the UTC at an inner radius of 45 cm and an outer radius of 84 cm, the 
RS consisted of both scintillation counters and embedded straw chambers, providing energy 
and range measurement of the charged particles, information on the ir + — > fi + — > e + decay 
sequence and measurement of photon activity. 
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1. Scintillation Counters 



The RS consisted of 19 layers of Bicron BC408 plastic scintillator, azimuthally segmented 
into 24 sectors as shown in Fig. [5J Layers 2-18 were 1.9 cm thick and 1.82 m long. Layer 
19 was 1.0 cm thick and was mainly used to veto longer range muons. Each scintillator 
in layers 2-19 was coupled through Lucite light guides to EMI 9954KB PMT's at both 
upstream and downstream ends. The innermost counters were 6.4 mm thick and 52 cm long 
trigger-counters (T-counters), defining the fiducial volume for charged K + decay products. 
The T-counters were thinner than layers 2-19 to suppress rate from photon conversions. 
Seventeen 1-mm-diameter WLS fibers (Bicron multi-clad BCF-92) with a pitch of 6.9 mm 
were embedded in each scintillator and coupled to a single Hamamatsu R1398 PMT at each 
end. 

The signal from each PMT of the RS scintillators was passively split 1:2:2 for ADC's, 
discriminators and fan-in modules, respectively. The discriminator output was sent to a TDC 
and the trigger. Each PMT was read by an ADC and a TDC. The TDC's (LeCroy 3377) 
recorded up to 16 hits in a 10.5 /is time window, and thus allowing for efficient detection 
of the fi + — > e + decay. The analog fan-in summed the signals from 4 PMT's on each end 



in the same hextant (four adjacent sectors) and 



single TD 



58| and was provided to mean-timers 



aver. This analog sum was read out by a 



63j for good timing on the photon veto in 



the trigger and a ^-measurement for each layer of the track. The mean-timer output from 
each layer in a hextant was ORed to provide input to the hextant photon veto algorithm in 
the trigger (vetoing more than one non-adjacent hextant). The TD's recorded the charge 
in 2 ns intervals (500 MHz sampling) in a 2.5 /j,s time window with a resolution of 8 bits. 
The 500 MHz sampling provided sufficient pulse shape information to separate pulses from 
different events as close as 5 ns apart, and enabled the detection of the ir + — > fi + decay as 
described in Section IIIIF3I The time window of the TD's was narrower than that of the 
TDC's in order to reduce the data size. 



2. Range Stack Straw Chambers 



Two range stack straw chambers (RSSC's) were located outside RS layer 10 and 14 60]. 



The inner (outer) RSSC consisted of two staggered layers of 24 (28) straws per sector with a 
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FIG. 6: Schematic end view of the inner RSSC. The tubes run through the beam direction. 

length 97.8 (113.0) cm. The average density of an RSSC was 0.054 g/cm 2 . Each straw tube 
was 3.4 mm in radius with a 50-/im-diameter gold-coated tungsten anode wire at the center. 
A schematic drawing is given in Fig. [6j The straw chambers were operated with 67% argon 
and 33% isobutane mixture with a trace of water in a self-quenching streamer mode at 3,450 
V. The local x-axis of a chamber was denned to be along the width of the chamber. There 
were a total of 48 chambers, with 2,496 straws installed in the E949 experiment. Because 
of access restriction, a pair of straws from the right and left halves in the same layer were 
connected at the upstream end to allow downstream-only readout. 

The position of the hit straws provided x — y position information, while the end-to-end 
time differences provided the z measurement. The z resolution of the E787 RSSC's was 
degraded due to a pulse- height dependent time- walk effect. For E949 an amplifier and two 
discriminators were installed on each channel to allow for low threshold timing discriminator 
and a high threshold logic discriminator (above the noise level). The z resolution was 
improved from 3.0 cm (RMS) observed in E787 to 1.5 cm (RMS). 

G. Photon Veto Counters 

The detection of activity coincident with the charged track was crucial for suppressing 
background processes that can mimic K + — > n + uu. Photons from K n2 and other radiative 
decays were detected by the hermetic photon system as shown in Fig. [H The photon 
detectors surrounding the K + decay vertex with a An solid angle coverage were located in 
the barrel, upstream and downstream end caps, and near the beam line. The photon system 
consisted of essentially every scintillator detector in experiment: the Barrel Veto (BV), the 
Barrel Veto Liner (BVL), the RS, the upstream and downstream End Caps (EC's), the 
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Upstream Photon Veto (UPV), the upstream and downstream Collar detectors (CO), the 
downstream Microcollar detector (/iCO), the Downstream Photon Veto (DPV), the RV, the 
IC, the VC, the AD and the target. The regions of the target, IC and RS traversed by the 
charged track were excluded from the photon veto. The AD and DPV were part of the E949 
detector upgrade but were only used in the 7ri/z/(2) analysis, where photon veto near the 
beam axis was more important. 

The 1.9- meter- long, 14.3 radiation length thick BV covering 2/3 of the 4tt sr solid angle 
was located in the outermost barrel region with an inner radius of 94.5 cm and an outer 
radius of 145 cm 5Jj . The BV was divided into 48 azimuthal sectors. Each sector consisted 
of four radial layers, in which there were 16 (innermost), 18, 20, 21 (outermost) layers of 
1-mm thick lead and 5-mm thick Bicron BC408 plastic scintillator. The light collected in 
the scintillators accounted for 30% of the total energy deposit in the BV. The azimuthal 
boundaries of each sector were tilted so that there were no projective cracks for photons 
from the decay vertex. Both ends of every module were read out by an EMI 9821KB PMT 
and the signals were recorded by an ADC and a TDC. The time resolution of individual BV 
counter was measured to be 1.2 ns. 

In order to improve the photon veto capability, the BVL, located between the RS and 
the BV, replaced the outermost layers 20 and 21 of the RS in E787. Each BVL counter was 
10 cm wide and 2.2 m long. There were 48 azimuthal sectors, each with 12 layers of 1-mm 
thick lead and 5-mm thick Bicron BC408 plastic scintillator, for a total thickness of 2.29 
radiation lengths. Both ends of the BVL modules were read out by EMI 9821KB PMT's 
and the signals were recorded by ADC's and TDC's. The eight adjacent sectors (hextant) in 
each end were grouped and read out by TD's. The timing resolution of an individual BVL 
counter was 0.7 ns. A comparison of radiation length coverage with and without the BVL 
is shown in FigJTl A factor of two improvement in the photon veto rejection of K n2 decays 
was expected from the BVL. 

The EC's had roughly one-third of the 4ir sr photon coverage 6JJ. The upstream EC 
detector consisted of seventy-five 25 cm long (13.5 radiation lengths) undoped Cesium Iodide 
(Csl) crystals segmented into four rings, and the downstream EC detector consisted of 68 
crystals in four rings. To maximize light collection the PMT's were directly coupled to 
the crystals through a Sylgard cookie formed over the PMT and a ultraviolet transmitting 
optical filter that selectively passed the fast component of the Csl scintillation light (with a 
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FIG. 7: Radiation length with (solid curve) and without (dotted curve) the BVL as a function of 
the cosine of the polar angle. These curves account for the contribution from RS, BV, EC and 
BVL. Other photon veto counters are not accounted for in this plot. 

decay time of a few tens of nanoseconds at a wavelength of 305 nm) and blocked the slow 
component. Since the PMT's were situated in the magnetic field, high-field PMT's were 
used 62] (Hamamatsu R5543 3" PMT's for the outer three rings and Hamamatsu R5545 
2" PMT's for the smaller inner ring crystals); the signals were split to ADC's, constant 
fraction discriminators (CFD), and CCD's. The CFD's outputs were ORed to provide an 
online photon veto signal and were also sent to TDC's. The pulses recorded in the CCD's 
were analyzed offline by a pulse-finding algorithm to provide the timing of the EC signals 
and to separate two pulses close in time (to reduce the possibility of accidental hits reducing 
the efficiency of photon detection). Since the EC was exposed to a high counting-rate 
environment near the beam line, good timing was required to reduce acceptance losses and 
the masking of i^photons from early accidental hits. 

The UPV mounted to the downstream face of the Cerenkov counter was 3.1 radiation 
lengths thick, with an outer dimension of 28.4 cm x 28.4 cm and an inner hole for the beam 
of 17.5 cm (horizontal) x 4 cm (vertical). It consisted of twelve layers of 2 mm thick Bicron 
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BC404 plastic scintillator and 1-2 mm thick lead plate. The scintillation light was read out 
by 12 x 21 WLS fibers coupled to two Hamamatsu R1924 PMT's. The UPV signals were 
sent to ADC's, TDC's and CCD's. 

The upstream (downstream) CO detector was located just upstream (downstream) of the 
EC's [ssj] . Both of the CO's consisted of twenty-four 2-mm thick lead sheets alternating with 
25 layers of 5-mm thick Bicron BC404 scintillator sheets, providing about 9 radiation lengths 
along the beam direction. Each scintillator layer was segmented into 12 wedges forming 12 
identical azimuthal sectors. Light from the wedges was readout by WLS multi-clad fibers 
(Bicron BCF99-29AA). Each wedge had 16 fibers glued in grooves in the BC404 scintillator 
layer. One end of the fiber was polished and aluminized to provide reflective mirror surface. 
All fibers for each sector (16x25=400) were terminated in a connector which was coupled 
to a 1.1 7-m long Lucite rod via a silicone cookie. The rods transmitted the light to EMI 
9954KB PMT's in a low field outside the magnet. Signals were sent to ADC's and TDC's. 
The CO's detected photons emitted with a small polar angle (0.970 < \cos6\ < 0.995) in 
both the upstream and downstream regions. 

The /iCO was installed between the inner face of the magnet end-plate and the target 
downstream of the downstream CO to give more photon detection coverage. The /iCO had 
an inner (outer) diameter of 15.6 cm (20.0 cm) and a length of 53 cm. The eight azimuthal 
sectors contained eight layers of 2 mm scintillating fibers (from the original E787 target 5l|) 
separated by seven layers of 60 /mi Pb. The 536 fibers from two adjacent sectors were read 
out by one EMI9954 PMT into ADC's and TDC's. 



H. Trigger 

The E949 trigger selected K + — > 7T + uu events from the large number of K + decays and 
scattered beam particles with requirements on the range of the 7r + track, the presence of a 
7r + — > p^v^ decay in the RS, the absence of other activity at the time of the ir + and the 
presence of a K + at an appropriately earlier time. 

The trigger was composed of two stages: a fast Level-0 trigger with a decision time of 
~100 ns and a slower Level-1 trigger with a decision time of 10-100 /as. The Level-0 trigger 
was based on signals from the beam, target, RS and photon veto systems, processed with a 
combination of commercial and custom-built boards. The Level-1. n trigger was composed of 
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two parts running in parallel, the Level-1.1 and Level-1.2 triggers, that involved the partial 
processing of TD and ADC data, and operated only on events that passed the Level-0 trigger. 

The original E787 trigger has been described previously 5l|. The trigger was upgraded 
for E949 with the addition of a new programmable Level-0 trigger board and mean-timers 
for the photon veto signals jf^ . 

E949 collected data with triggers for the K + — > n + uu signal and for other physics studies 
and calibrations. In this paper, the triggers for calibration and other physics are referred to 
as the "monitor" samples. 



1. Trigger Architecture 

For physics and most calibration triggers, Level-0 required a signal from the beam system 
and a charged track in the RS. A beam K + or ir + was identified in NIM logic by a coincidence 
of hits from the Cerenkov counter, B4 (BA), target (Target) and AGS spill gate (Spill). 
This KB signal served as the Beam Strobe (BS) for the trigger. A charged track T • 2 was 
identified by a coincidence of the IC and T-counter and second layer (mean-time of the two 
ends) in a single RS sector. Additional requirements were placed on the crude estimate of 
the track range, which was the track length derived from the deepest layer of the RS hit in 
coincidence with the T • 2 from that sector or the two proceeding sectors (for a positively 
charged track). Another requirements were placed on the photon veto (EC, RS, BV, BVL) 
and the delayed-coincidence (DC). The DC was an OR over the six coincidences of the 
individual IC's and the delayed KB signal such that a track in the IC must typically have 
occurred at least 1.5 ns later than a prompt signal. The T • 2 signal, an OR of the 24 
individual sector T ■ 2's, typically served as the Detector Strobe (DS) that gated all of the 
ADC's not associated with the beam system or target and provided a common stop for many 
of the TDC's; the T • 2 signal introduced 40 ns of dead time. When an event passed the 
Level-0 trigger, the dead time was extended to 100 ns to allow further processing. Additional 
triggers that did not contain both a beam signal and a T • 2 were produced on a separate 
trigger board that fed directly into the trigger bus (with external coordination of the dead 
time). 

The Level-1.1 trigger (Ll.l) was based on information from the TD system for the hex- 
tant containing the RS counter in which the charged track was determined to have stopped. 
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Stopped 7r + 's were preferentially selected (over /i + 's) by looking for the n — > n decay by 
comparing the pulse height to the pulse area (for early decays) or by looking for a second 
detached pulse. This was done by a custom-built ASIC which had access to the TD memo- 
ries. The ASIC could also reduce the TD data size for readout by discarding the waveform 
samples outside the "prompt" time window, keeping instead a calculation of the pulse area 
and leading-edge time. The "prompt" window typically extended from 0.5 (2) /is before 
(after) the ir + track. Level-1.1 typically provided a decision in about 10-20 /is. 

The Level-1.2 (LI. 2) used data digitized by the ADC's and had two components. One of 
them, the "Level-1.1 afterburner" rejected events with an accidental hit near the stopping 
counter. Such hits might defeat the Level-1.1 with an apparent double pulse. The other 
component, the ll HEX afterburner", was used for the photon veto, rejecting events with 
hits in both of the two adjacent hextants when the T • 2 counter and the stopping counter 
were found in the same hextant. This usually indicated that a single hextant fully contained 
the trajectory, thus those observed in the other hextant came from another particle. The 
Level-1.2 trigger introduced a dead time of up to 100 /is per Level-1.1 trigger. 

2. K + — ► ■k + uv Triggers 

The triggers as described below were valid for the bulk of the 2002 running period. The 
trigger conditions for K + — > 7i + uu were designed differently according to the 7r + momentum. 
For high ir + momentum, the irvv(l) trigger condition was defined as: 

Tzvu(l) = KB ■ DC ■ (T -2) ■ (Q ct + 7 ct ) ■ (15) 
19^ • zfrf ■ LOrrl ■ (BV + BVL + EC) ■ 
HEX ■ Ll.l ■ Ll.2, 

where 

KB = C K ■ BA ■ Target ■ Spill. (16) 

The 6 c t + let required that the ir + reached the 6 th or 7 th layer of the RS and suppressed the 
copious 3-body K + — > 7r + 7r~7r + and K + — > 7r + 7r 7r° backgrounds. The 19 c t signal required 
that the charged track should not reach the 19 th layer in order to suppress K^ background. 
The "ct" designated the RS sectors that were associated with a T-2 (T-2 sector plus the next 
two clockwise sectors: this was the direction that a positive particle moved in the magnetic 
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field). The zfrf condition was a fiducial cut on the ^-position of the charged track in each 
layer, vetoing tracks that exited the fiducial volume. The LOrrl was a refined calculation of 
the charged track range, which included the number of target fibers hit and a measurement 
of the z-position of the track from flash TDC's on layers 3 and 11-13 as well as the deepest 
layer of penetration of the track; this rejected events with long range such as the /i + from 
Kya. decay. The photon veto BV + BVL + EC and HEX were from the BV, BVL, EC and 
RS, respectively, which removed events with photons such as K n2 , and K^y 

The number of K +, s which met the KB trigger requirement when the detector was live 
(Nk), was recorded at the end of each AGS spill. The total exposure from the ttvv(1) trigger 
data stream was 

N K = 1-77 x 10 12 . (17) 

E949 also provided a t[uv{2) trigger for lower momentum 7r +, s. Since the study involves 
different background mechanisms and a different kinematic region, the result will be pre- 
sented in a separate paper. 



3. Monitor Triggers 

In addition to the K + — > tt + vv triggers, there were monitor triggers for calibration and 
normalization, as well as triggers for other physics modes. All triggers were prescaled to 
reduce the impact on the total dead time except for the K + — > n + uu triggers and the "7" 



trigger which required the 



and K + — > 7r + 7 decays 66]. 



Dresence of photons in the barrel region for studying K + — > 7r + 77 



To monitor detector performance several processes were employed, including K n2 and 
decays, beam particles scattered into the detector fiducial volume, beam K + , charge 
exchange events and cosmic rays. These monitor samples were used for calibration and 
acceptance studies. They were taken simultaneously with the ivvv(l) and itvv{2) triggers 
to reflect any condition changed into the acceptance and background calculations. 

The Kyu decay has the largest branching ratio. Since the final state does not contain 
photons or additional tracks and the daughter fi + does not interact strongly, it was a con- 
venient sample for a variety of acceptance measurements as well as several calibrations. 
This sample was also used for the normalization of the experiment; the measurement of the 
K^2 branching ratio effectively normalized the counting of K + stops to the well-known 
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branching ratio. The K\i2 trigger condition for K^ was defined as follows 

A>2 — KB • (T ■ 2) • (6 d + 7 ct ) ■ (17 ct + 18 ct + 19 ct ). (18) 

The final state of the K n 2 decay mode contains one charged 7r + and two photons from 7r° 
decay. Since the tt + momentum is monochromatic, the K n2 sample was used to check the 
measurement of the charged track momentum, range and energy. Also the 7r + 's were used 
to study particle identification, while the photons were used to study the photon veto. Two 
K n 2 triggers were defined, a loose one, Ktt2(1): 

Kn2(l) = KB • (T • 2) • (6 d + 7 ct ) ■ 19^. (19) 

and a tighter one, A7r2(2): 

Att2(2) = Avr2(l) • HEX ■ Ll.l ■ L1.2. (20) 

The Ktt2{2) monitor trigger was used only for calibration in this analysis. 

Among the incoming beam particles there were many 7r +, s, including some that scattered 
into the fiducial volume of the RS. These 7r scat events were identified as 7r + 's by the C n and 
had an in-time track in the RS. The kinematic features of this 7r + sample were almost the 
same as the K + — > tt + uu signal except that the target pattern was different. This sample 
was suitable for calibrating the ionization energy loss (dE/dx) of 7r + 's and for studying the 
acceptance. The trigger condition for scattered beam particles was defined as: 

7i scat = nB ■ DC ■ (T • 2) ■ (6 ct + 7 a) ■ (BV + BVL + EC) ■ HEX. (21) 

where 

nB = C 7T -BA- Target ■ Spill. (22) 

The charge exchange process (CEX), K + n — > p+K° followed by K\ — * tx + 1^vi can mimic 
K + — > 7i + uu events when the charged lepton l~ from the A° decay had a low momentum 
and was undetected. The largest uncertainty in this background was the determination of 
the reaction rate. Since K° has an equal fraction of A° and K§, the K + n — > p + K$ process 
with a K^, decay to n + 7T~ can be used to measure the reaction rate for the determination of 
the background from K + n — > p + K\. Since the Kg — > 7r + 7r _ could be cleanly identified, a 
CEX trigger with two charged tracks was defined as: 

CEX = KB-DC- 2(T • 2) • (6 d + 7 ct ) ■ EC + nB. (23) 
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Type 


Model 


Standard 


Resolution 


Subsystems 


ADC 


LRS 4300B 
LRS 1881 


CAMAC 
Fastbus 


10 bits 
13 bits 


RS,BV,BVL,EC,Beam 
Target,UTC 




LRS 3377 


CAMAC 


0.5 ns 


RS, BVL 


TDC 


LRS 1879 


Fastbus 


2 ns 


UTC,BV,Target 




LRS 1876 


Fastbus 


1 ns 


EC,RSSC,Beam 


WFD 


TD 

CCD 


Fastbus 

I 

Fastbus 


500 MHz sampling 
i bits, up to 10 /j,s depth 
500 MHz sampling 
8 bits, 256 ns depth 


RS,BVL,IC 
Beam, Target, EC 



TABLE I: Digitizing electronics for E949. 



In addition, we also denned triggers for beam K + and cosmic ray, which were used in 
trigger efficiency measurement and detector geometrical alignment. 



I. Data Acquisition 

Analog- and discriminated- signals from the detector were digitized by commercial ADC 
and TDC, and custom-built waveform digitizer (TD and CCD) systems. When an event 
was accepted by the trigger system, the digitized data for the event were transferred to a 
buffer module or a local crate controller. At the end of each spill, the data for the spill were 
transferred to a host computer. A summary of the digitizing electronics is shown in Table [H 



For the Fastbus systems, SLAC Scanner Processor (SSP) modules (67J served as crate 
controllers and also to read out, reformat and buffer the data from the front-end after each 
trigger accepted. The CAMAC ADC's were read out through the FERA bus by a Struck 370 
QDP DSPfFastbus) module. The CAMAC TDC's were read out by custom-built DYC3 
modules 68| which pushed the data into VME memory boards. The readout time per event 
(as determined by the slowest crate) was typically 850 fis. 

At the end of each spill, the data from the Fastbus buffer memories were read out via the 
cable segment (12-15 MB /sec) by Struck 340 SFI modules, each controlled by a MVME 2604 
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single-board computer (SBC) running VxWorks. The VME memory boards were read out 
by a separate SBC. Data were transferred from the SBC's to the host computer (SGI Origin 
200) via Ethernet (9 MB/sec per link) through a simple network switch. Event fragments 
from the readout segments were combined by Event Builder processes running on the host 
computer. Complete events were distributed to "consumer" processes which included data 
logging and online monitoring. The K + — > tt + uu triggers were written to two DLT-7000 
drives at 5 MB/sec per drive; a third DLT drive was used to log monitor triggers. 

A slow control system, based on the MIDAS 69| framework, ran independently of the 
main DAQ system and was used to monitor a variety of experiment conditions, including 
crate voltages and temperatures. 

Under typical running conditions, we wrote ~300 events per spill with a typical event 
size of ~80 kB. This was well within the maximum throughput of the system of about 
50 MB/spill. The DAQ dead time was due entirely to the speed of the event-by-event 
readout of the front-end electronics at the crate level. The total dead time introduced by 
the trigger and DAQ was typically 26%. The E949 experiment collected its physics data 
for 12 weeks from March through June of 2002 or about 20% out of the total beam time 
approved by DOE's Office of High Energy Physics. 



III. DATA ANALYSIS 

The branching ratio of the K + — > ir + uu decay in the SM is ~10 -10 as discussed in Sec- 
tion [H Unlike the and K n2 backgrounds, the K + — > tt + vv signal is continuous with no 
peak. To establish that any possible observed candidate event was really from K + — > ir + uu, 
we required that the backgrounds were suppressed to a level substantially below one event 
with small uncertainty. All the detectors were calibrated before the background and accep- 
tance studies using Kfi2, Ktt2(1), Ktt2(2) and cosmic-ray and beam trigger events. After 
a brief overview of the background sources, the data analysis technique will be described: 
including selection criteria, background evaluation, acceptance measurement and the signal 
candidate search. 
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FIG. 8: Range in plastic scintillator (cm) versus the momentum (MeV/c) of the charged particles 
for events that passed the irvvil) trigger. These data represent ~0.3% of Nk- 

A. Overview of Background 

Data selected by the irvvil) trigger were primarily from background events as shown 
in Fig. [HI These events were classified into stopped- _fT + -decay-related and beam-related 
backgrounds. 
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1 . Origins of Stopped-K + -Decay Background 



The stopped- -?T + -decay backgrounds were categorized into two types: 7r + -related, 
related backgrounds. As can be seen in Fig. |U multi-body K + decays were suppressed by 
setting a signal momentum region higher than the K n2 peak but lower than the peak. 
However, since K^i and have such large branching ratios (20.92% and 63.44% 181]). 
migration into the signal region through either resolution or scattering effects was a signif- 
icant background. Radiative decay and K + — > n + it°v {K^) decays accounted for the 
majority of the observed fi + band events in Fig. [HJ For the decay mode, the background 
originated from a peak event or a K^ 2 low momentum ("tail") event if the particle 
identification was fooled. This was also applicable to the K n2 decay mode when the pho- 
tons escaped detection. Because of phase space limits, the 7rvv(l) analysis only needed to 
consider the _fT + -decay-related backgrounds from the K^, H + band and K n2 decay. 



2. Origins of Beam Background 

The beam-related backgrounds were categorized into three types: single-beam back- 
ground, double-beam background and CEX background. The first two beam backgrounds 
accounted for most of the ir + band shown in Fig. [8j 

The single-beam background consisted of the following components. (1) A K + entered 
the target and decayed in flight to a ir + plus a tt° as illustrated in Fig. The kinematic 
values of the ir + were shifted upward to the signal region by the Lorentz boost, faking a 
K + — > ir + uu signal event. (2) A ir + in the beam was misidentified as a K + , scattered in 
the target and entered the fiducial region of the detector as illustrated in Fig. [9j This 
7r + (referred to as scattered n + ) could mimic the target fiber pattern for signal and have 
kinematic values in the signal region. Rejection of these two background types required both 
good 7r + / K + beam identification and time resolution for delayed coincidence measurements. 



The following cases were classified as the double-beam background. A K + came to rest 
in the target accompanied by another K + entering the target that decayed in flight to a tt + , 
which traversed the fiducial region of the detector (Fig. [TOl) . The second case was similar, 
except that the beam 7r + scattered in the target and entered the fiducial region (Fig. [TOl) . 
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FIG. 9: Schematic diagrams of the single beam background: 1) single beam K + background and 
2) single beam tt + background. The various detector elements and acronyms are described in the 
text. 
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FIG. 10: Schematic diagrams of the double beam background: 1) double beam K + — K + back- 
ground and 2) double beam K + — tt + background. The various detector elements and acronyms 
are described in the text. 

Both cases could imitate a K + — > tx^vv signal if the decay products from the first K + were 
missed. Rejection of these two backgrounds relied on the ability to observe extra activity 
the beam instrumentation, target and RS coincident with the delayed decay. 
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FIG. 11: Schematic diagram of the charge exchange interaction background. The various detector 
elements and acronyms are described in detail in the text. 



The CEX background could occur if a K + produced a K° in the target and if the K° 
turned into a K\ that subsequently underwent semileptonic decay. This process could 
produce background if the charged lepton from a K\ decay was undetected and the ir + 
satisfied the kinematics of the signal region (Fig. [TT]) . Rejection of the CEX background 
was achieved by using the fact that a K\ usually did not deposit energy along the path in 
the target, leaving a gap observed between K + and 7r + fibers in the target. Also exploited 
was matching between the reconstructed ^-position of the decay and the energy deposited 
by the incoming K + and the short flight time of K\ in the target. 



B. Analysis Method and Strategy 

Disentangling the K + — > tt + vv decay from background in this experiment was challenging 
due to the poor kinematic signature and very small expected rate of the K + — > tt + uu signal. 
These necessitate enormous suppression of background events by rejecting events with very 
low levels of extraneous activity. This high level of veto makes it impractical to accurately 
simulate the background rejection power of the detector at the required ~ 10~ 10 level of 
sensitivity. Therefore, an accurate estimate of the background must be obtained from the 
data. 
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1. Blind Analysis Method 

A "blind" analysis method was developed to search for the K + — > tt + uu signal in data 
samples. In this method, background sources were identified a priori. The signal region 
for the irui>(l), determined so that the sensitivity was optimized, was "blinded" or hidden 
until the background and acceptance analysis was completed. When possible, selection 
criteria were developed using the monitor samples to avoid examining the signal region. If 
monitor samples were inadequate and the 7ri/z/(l) trigger sample were required, at least one 
selection criterion distinguishing signal from backgrounds was inverted (i.e. used to select a 
background region) to avoid examining the signal region. In addition, the final background 
estimates were obtained from different samples than that used to determine the selection 
criteria. Each set of three consecutive events were selected for "1/3" and "2/3" 

sample groups. The 1/3 group was used to determine the selection criteria and an unbiased 
background estimate was obtained from the 2/3 group. The signal region was examined 
only after the background analysis was completed. 

2. Bifurcation Method for Evaluating Background 

The principal method for background evaluation relied on information from outside the 
signal region and involved the application of two complementary but uncorrelated cuts. 
Fig- E±] illustrates this bifurcation method showing the parameter space of two cuts, "CUT1" 
and "CUT2". The number of background events in the signal region (i.e., region "A") was 
A events. If the two cuts are uncorrelated, that is, if the rejection of a cut does not depend 
on the rejection of the other cut, the ratio of the number of background events in region 
"A" to region "B" must be equal to the ratio in region "C" to region "D", i.e., A/B = C/D. 
Background events in the signal region are therefore obtained from the relation A = BC/D. 

In practice, the present bifurcation analysis was done through two branches. A "nor- 
malization branch" analysis was used to obtain the number of events in "B" region. A 
"rejection branch" study was done to get the ratio of D/C. The rejection was defined 
as R = (C + D)/C. The background level in the signal region was then estimated as 
Nskg = B/(R — 1). For the case of very small statistics, "CUT1" was subdivided into two 
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FIG. 12: Pictorial explanation of the bifurcation method. Background level in Region "A" can 
be estimated from the observed number of events in the other regions, if CUT1 and CUT2 are 
uncorrelated. 

cut categories, and B was estimated in the same way as the "first" bifurcation. 

To check if the two bifurcated cuts were uncorrelated, cuts were loosened simultaneously 
(as described in Section IV A). The loosening factors were controlled by the predicted back- 
ground functions, in which the loosening factors were the inputs and the outputs were the cut 
positions. By design the functions provided the relative background level and the acceptance. 
The background level provided by the functions should agree with the observed number of 
events within the newly defined regions if CUT1 and CUT2 were uncorrelated. This method 
thus provided input to the evaluation of systematic uncertainties (Section III! G 101) . 

3. Analysis Strategy 

The data analysis used the following key steps to determine the selection criteria, evaluate 
the background level, investigate the systematic uncertainty, measure the acceptance and 
finally obtain the branching ratio. 
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Data were first reconstructed and processed with a number of selection criteria to 
remove obviously bad events. Then the surviving events were divided into 1/3 and 
2/3 portions, in which three sub-samples were also skimmed out according to different 
background features (Section IIIIEf) . 

Blind analysis was adopted in designing, calibrating and tuning the selection criteria. 
Signal-like and background samples were taken from the monitor trigger samples when 
applicable. When the three skimmed sub-samples were used to have the same exper- 
imental features of the signal, at least one critical selection criterion was inverted to 
avoid examining the signal region (Section IIIIFj) . 

The background level was initially evaluated by applying the bifurcation method to 
the 1/3 data and controlled to be much less than one event by tuning the selection 
criteria. At least two uncorrelated cuts with large background rejections were chosen, 
data samples (Section IIII G llUlI G 5j) . 

The final cut positions were optimized with respect to the signal to background level 
estimated from the 1/3 data. This was achieved using the predicted background 
functions estimated in the bifurcation analysis (Section IIII G 6tHlI G 7j) . 

Correlations between the bifurcated cuts were checked by conducting a series of re- 
evaluations of the background levels outside the signal region. A study of single cut 
failure was also conducted to investigate a possible flaw in the technique. To avoid 
the potential bias of the 1/3 background study, the final background evaluations were 
taken from the 2/3 portion fSection lITrG"8lllllG9p . 

Acceptances were measured with the monitor trigger samples except for the signal 
phase space, the trigger efficiency and those that could not be extracted from the 
monitor trigger data. These were obtained from Monte Carlo. The branching ratio 
of K W 2 was used to validate the acceptance measurement. Single event sensitivity was 
obtained from the total K + exposure and the acceptance. (Section IIIIHj) . 

The signal region was examined. Events observed in this region were all considered as 
the signal candidates (Section IHIIj) . 
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• The branching ratio was obtained from a likelihood analysis incorporating the pre- 
dicted background rate and acceptance within the signal region (Section HVl) . 

C. Track Reconstruction 

Throughout this analysis, the events were reconstructed under the assumption that they 
were K + — > tt + i>v events with only a single 7r + track in the detector. 

1. Beam Time Measurements 

To fully reconstruct an event, the initial time of the beam particle was required. The 
beam instruments provided several beam times from TDC and CCD measurements on the 
Ck-i the CV, the BWPC and the B4. For the TDC measurements, offline analysis treated 
all PMT hits coincident with each other as a cluster. The average time of the TDC hits in 
each cluster gave tc K , t^, tsw an d tBM, respectively. The CCD measurements were used 
to discriminate cases with more than one particle in a beam (referred to as pileup). 

2. Clustering in the RS 

The track reconstruction routine started by finding clusters in the RS. The hit counters 
of a positively-charged track (the track counters) were searched for using the TDC timing 
information in the RS. A good T ■ 2 sector was found as that closest in time to the DS. From 
this T ■ 2 sector, adjacent counters in subsequent layers within 10 ns of the T • 2 time and 
with energy greater than 0.5 MeV were selected in the outgoing and clockwise direction to 
find the track counters. A good cluster involved at least six consecutive layers from inside 
to outside. Once an RS track cluster was identified, the track time (t rs ) was computed 
by averaging all the time measurements of the track counters. The stopping counter was 
defined as the one in the outermost layer which was in the most clockwise direction. The 
T • 2 sector in a cluster served as a guide for tracking in UTC. 
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3. Tracking in the UTC 



When an RS cluster was established, the UTC tracking routine started searching for 



clusters of hit anode wires in the x — y plane [59[. Hit wires in each super-layer were 
grouped into clusters based on their spatial proximity to one another. A straight line fit 
provided a crude vector in each super-layer. The solutions due to the left-right ambiguity 
were included at this stage. These vectors were then linked to form a track segment in 
the UTC. A circle fit was performed with a set of drift distances with left-right ambiguity 
resolution. The radius of the circle gave the measurement of transverse momentum. 

If a UTC track in the x — y plane was found, the corresponding track projection on the 
cf) — z plane (<f) was defined as the revolution angle with respect to the closet approach to 
the vertex) was then sought in the clusters of hit strips on the UTC cathode foils. A time 
window of ±15 ns was used to reduce the accidental hits. The calculation of the z position 
for a cluster adopted the ratio method suggested by Ref. [3], which used the three strips 
with the highest ADC counts to derive the centroid and to reduce the bias on the z position 
measurement. A straight line fit was performed in — z plane if the z hits were found in at 
least 3 foils, thus determining the slope and intercept. The slope was then used to convert 
the measured transverse momentum to the total momentum. 

In case of more than one track pointing to the same T • 2 sector, a good UTC track was 
defined as the one closest to the first RS sector crossing point or to the clock-wise edge of 
the stopping counter otherwise. The UTC tracking efficiency was checked with the Kfi2 
monitor data and was measured to be better than 95%. The inefficiency was due to the fact 
that for the high rate environment accidental hits in the UTC caused problems in pattern 
recognition. 



4- Tracking in the Target and T>4 

Traveling almost parallel to the target fibers, the incident K + deposited significant energy 
in each fiber (usually > 4 MeV) and was in coincidence with the BS. The daughter 7r + 
traveled nearly perpendicular to the target fibers, and thus left less energy in each fiber 
(1 MeV on average) and was in coincidence with the t rs . The BS, t rs and fiber times, 
positions and energies were the key elements to identify the K + and 7r + fibers. All the K + 
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and 7r + fibers were linked to form a K + cluster and a ir + cluster. A good event should only 
have one K + cluster and one ir + cluster. 

After a UTC track was found, the target pattern recognition routine started to look for 
the fibers belonging to the K + path and the ir + path separately on a 1 cm wide strip along 
the UTC track extrapolated into the target. The corresponding energies for K + and 7r + [Ex 
and Etg) and times (tx and t n ) were calculated from the sum and average in the clusters, 
respectively. The range of ir + in the target (Rt g ) was calculated as the helix traversed by 
the 7r + from the K + decay position to the inner surface of the IC with the polar angle 6 
correction. 

Because of ambiguity in the entry and stopping ends in the pattern of the incoming beam, 
the target reconstruction routine used the B4 measurement on the K + entrance point, which 
was reconstructed by clustering the hits in the two B4 layers to determine the beam time 
t BM , the energy- weighted beam position and the energy loss in B4. The cluster with t BM 
closest to the tx was chosen as the one caused by the K + beam. A 0.36 cm position 
precision was obtained by the B4 hodoscope in the x — y plane. With the B4 position 
measurement, the target reconstruction was repeated to give a better pattern recognition. 
The K + decay vertex in the x — y plane was determined by the K + fiber closest to the UTC 
track but furthest from the K + entrance point determined by the B4, while the z position 
was calculated from the UTC track extrapolation in the <p — z plane. 

Target CCD information improved the pattern recognition and ir + energy measurement 
when pileup occurred in a K + fiber. When tx and t n were separated by more than 2 ns, the 
CCD pulses in all of the K + fibers were studied to identify any hidden ir + energy. 

Isolated hit fibers outside the 1 cm search strip and within the time window of t rs ± 5 ns 
caused by possible photon(s) or other beam particles were classified into "photon" fibers if 
their energies were greater than 0.1 MeV. 

Once the first round of target reconstruction was finished, the ir + passage and K + decay 
vertex were determined in the target. The procedure of UTC track fitting was repeated 
with this additional information. This aided the resolution of the left-right ambiguity in the 
UTC track reconstruction. Iteration of the target reconstruction was also performed with 
the improved UTC track. 
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5. Track Passage in the IC 



An allowed IC hit pattern was either one or two adjacent sectors per event. The IC 
provided energy loss (Ejc) and time (tic) measurements. If there was an IC sector crossing, 
the Eic was from the sum of two IC sectors and the tic was from the energy-weighted 
average. It was observed that 1% of tracks had extra hits in the IC's, confounding the 
measurements of energy and time in IC's. This was resolved by using TD information in 
addition to the TDC and ADC information. The tic was always the one closest to the 
t rs and, the corresponding energy was taken as Ejc- The range Ric was computed as the 
length of the extrapolated UTC track from the inner to the outer IC radius. 

6. Tracking in the RS and RSSC 

With a charged track reconstructed in UTC and target, the tracking in RS started from 
the previously found RS cluster. The stopping counter was first analyzed by fitting the TD 
information with a double pulse assumption to find a n + — > fi + decay signature. This also 
determined the /x + energy (E^) from the tt + decay at rest. 

A sector crossing point as illustrated in Fig. [13] was searched for in the RS. A 
K + — > it + uu candidate should not have more than 2 sector crossings. Precise position mea- 
surement in x — y plane was obtained from the sector crossing points. The z positions were 
determined by using the end-to-end time differences in the hit counters except the T-counter. 
The average z position resolution was observed to be about 4 to 5 cm. 

Another precise position measurement was provided by the RSSC's. All the adjacent hit 
straw chambers in the RSSC of the hit sector were grouped to form a cluster. The candidate 
cluster in each sub-layer of RSSC was defined as the one closest to a series of arcs drawn 
from the UTC track extrapolation to the T counter through the sector crossing point (s). 
The x — y position was obtained from the average of two sub-layers in the same RSSC. To 
minimize the effect of cross-talk, the earliest hit in one sub-layer was chosen as the true hit. 
The z measurement was given by the time difference, and the precision was about 1.5 cm. 

An RS track fit in the x — y plane used the entrance point provided by the UTC track 
extrapolation, the sector crossing point (s), the RSSC hit position(s) and the expected path 
length predicted by the energy losses in the RS layers, taking into account 7r + track propa- 
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FIG. 13: Illustration of the Range Stack track fitting in x — y plane. All the relevant definitions 
are given in this plot. The arc represents the extrapolation of the fitted UTC track. 

gation with energy loss given by the Bethe-Bloch equation in a 1 Tesla magnetic field. The 
X 2 of the fitted track was minimized by changing the incident momentum and the angle at 
the entrance to the RS (Fig. [T3~j) . 

The total energy loss of the track in the RS (E rs ) was obtained by summing up all the 
energy losses in the track counters, with from the fit to pulse shape subtracted from the 
energy in the stopping counter. The range in the RS (R rs ) was calculated from the path 
length of the fitted track with the polar angle correction. The range in the stopping counter 
was estimated from the 7r + energy loss. 



7. Kinematic Measurements of a Track 

The total range R and energy E of the track were calculated as R = R tg + Ric + Rrs and 
E = E tg + Ejc + E rs , respectively. The total momentum P of the charged track was obtained 
from the UTC with corrections due to energy loss in the target and IC. Since the momentum 
reduction in both the target and IC were calculated from the range, some correlation between 
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P and R was expected, and thus they should not be treated as uncorrelated in the bifurcation 
analysis. All of these three kinematic measurements also included tiny contributions from the 
inactive materials in the UTC. In this analysis, the momentum, energy and range resolutions 
were measured to be 1.1%, 2.8% and 2.9% (RMS), respectively, from a study of K n2 events 
(Table OH). 



D. Monte Carlo Simulation 



Detector responses were modeled by a Monte Carlo simulation package, which was de- 
veloped for the E787 experiment and maintained or improved for E949. The pa ckag e used 
several subroutines from the electromagnetic-shower simulation package EGS4 71| and a 
number of routines written specially for the experiment, including all of the detector ele- 
ments, except for the beam instrumentation upstream of the target. Simulation samples 
were generated with the same format as the data except for omission of the pulse-shape 
information and most of the beam counter information. 



1. Simulation of K + Propagation 

The simulation of K + propagation started from a beam file, which contained K + events 
with a list of measured parameters: the K + stopping position, the tx, the number of K + 
fiber hits, the number of accidental fiber hits, the B4 hit position, the stopping target fiber 
element, the time, energy, fiber element for each K + and accidental hit in the target. This 
file was obtained from an analysis of the stopping distribution of the Kfi2 monitor data. 
With this beam file and the corrections described in Section HHP 31 the K + propagation had 
exactly the same target patten of the data in the simulation. Every K + decay started from 
the stopping fiber at the given K + stopping position. The K[i2 monitor data had sufficient 
statistics for all the simulation studies in this analysis. 



2. Simulation of K + Decay Product 

The K + — > tx + uu decay was generated with the matrix element of semileptonic 
decay via the V-A interaction, while the K^ 2 decay and the K v2 decay were generated via 
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pure phase space. Among the K + decay products, photon and electron interactions and 
their energy deposits were calculated using the routines from EGS4. For charged particles, 
the energy deposits were calculated by adding the energy losses of each ionization along the 
steps taken by the particles. The number of ionization and excitation events was determined 
by dividing the total average energy deposited along the step, obtained using the Bethe- 
Bloch formula, by the minimum energy that a particle lost in a collision. Multiple Coulomb 
scatterings of charged particles with various nuclei in the detector were calculated according 



to the theory of Moliere 
form factor of the nucleus 



with corrections for the spin of the scattered particle and the 



731 ] . Hadronic interactions of positively charged 7r + 's in the plastic 
scintillators were calculated using a combination of data and phenomenological models [7^] . 
An option in the simulation package allowed users to turn off nuclear absorption reactions 
and decays in flight of 7r +, s. This was useful in the study of the acceptance. 



3. Simulation of Trigger 

All the trigger conditions were simulated except for the DC, Ll.l, LI. 2 and those related 
to the beam instruments. Since the beam file generated from the K^2 monitor data included 
the trigger bias, a correction was needed to recover the true beam distribution. This was done 
by introducing a weight function derived from a comparison between simulated K\x1 triggers 
and data. In order to get a better precision, perpendicular tracks were used to calculate the 
K\i2 trigger acceptance as a function of radial distance of the stopping position in the x — y 
plane. 



4- Comparison between Data and Simulation 

The performance of the Monte Carlo simulation was checked by comparing kinematic 
resolutions between data and simulation as given in Table [Til for the K v2 decay mode. A 
0.15 cm deviation was observed in the range resolution, which is still not understood. This 
difference could have affected the acceptance estimate and was investigated when performing 
the acceptance study (see Section lHIH5l) . It should also be emphasized that the main role of 
Monte Carlo simulation in the E949 experiment was to estimate the acceptance factors that 
could not be obtained from real data (e.g. geometrical and relevant trigger acceptances). 
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Experiment ap (MeV/c) (Jr (cm) oe (MeV) 
Data 2.299± 0.006 0.866±0.002 2.976±0.005 

MC 2.399± 0.029 1.018±0.008 3.018±0.025 

TABLE II: Comparison of the momentum, range and energy resolutions between data and simu- 
lation for the K n 2 peak. 

E. Data Processing and Pre-selection 

Data were stored on the DLT tapes with a total data size of 7 Tera bytes. Two steps 
of data processing ("Pass 1" and "Pass 2") were taken to reduce and skim the data to a 
reasonable volume. 

1. Pass 1 

Pass 1 involved filtering cuts, which consisted of event reconstruction quality cuts and 
loose /i + background rejection cuts. Runs with trigger or hardware problems that could 
not be corrected offline were removed from the data analysis. Tracks were required to be 
successfully reconstructed, not to stop in those detector elements which had known hardware 
problems and to have momentum < 280 MeV/c. 

A 7r + — > fi + double pulse was required to be found by the fit to the TD pulse in the 
stopping counter and no extra hits were found in the other 3 sectors associated with the 
stopping counter TD channel. Since photon activity around the stopping counter could cause 
confusion in the energy measurement, events with sector crossing in the stopping layer were 
rejected. Also rejected were the events with a charged track that came to rest in the support 
materials for the second RSSC layer embedded between the 14 th and 15 th RS layers. Pass 1 
reduced the data volume by a factor of two. 

2. Pass 2 

Pass 2 involved cuts that were applied to the sample of events surviving Pass 1 to skim 
the data into three categories according to the event features. Each category was uniformly 
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divided into a 1/3 portion and 2/3 portion. Pass 2 consisted of five data skimming criteria: 
quality of the target reconstruction, loose photon veto, quality of the 7r + — > fi + double 
pulse fitting result, single beam K + requirement and delayed coincidence cut. These are as 
discussed below and listed in Table IIHI 

• The target reconstruction required that the K + decay vertex be inside the target 
volume, \tx — £bm| < 4 ns, \t n — tic\ < 5 ns and Eic be consistent with that expected 
from the calculated n + range. 

• The loose photon veto rejected events for \t — t rs \ < 2 ns with an energy in the BV 
greater than 1.5 MeV, \t — t rs \ < 1.5 ns with an energy in the EC greater than 3.5 MeV, 
\t — t rs \ < 1.5 ns with Ers > 3.0 MeV, or \t — t rs \ < 1 ns with E tg > 5.0 MeV, where 
t was the time measurement in each photon veto counter. 

• The 7r + — > yU + sequence required a /i + decay pulse in the stopping counter and the 
absence of hits within ±2.5 ns of the /i + time around the stopping counter. 

• The beam requirements for signal events were such that the energy loss in the B4 was 
greater than 1 MeV, the tsm differed from t rs by more than 1.5 ns, and the number 
of C n hits was less than 4 with \tc v — t rs \ < 1.5 ns. 

• The delayed coincidence required t n — tx > 1 ns. 

As defined in Table [TTT| Skim 1(4), 2(5) and 3(6) were enhanced in K^, fi + and beam 
backgrounds, respectively, for the 2/3(1/3) portions. These six skims of Pass 2 output data 
facilitated the studies to develop the final selection criteria and evaluate the backgrounds. 
Signal candidates, if any, remained in the six streams. Study of selection criteria using these 
samples was always done with at least one cut inverted to ensure that a blind analysis was 
conducted. 

F. Selection Criteria of post Pass 1 and Pass 2 

Further selection criteria were designed and applied in order to gain more background 
reduction. According to their characteristics, the selection criteria were classified into four 
categories: single beam K + selection criteria, kinematic reconstruction, ir + identification 
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Cuts 



Skim 1 & 4 Skim 2 & 5 Skim 3 & 6 



Target reconstruction 
Photon veto 
7r + — ► fi + sequence 
Beam K + requirement 
Delayed coincidence 

TABLE III: The cuts used in Pass 2 for selecting the data streams. Skim 1(4) , 2(5) and 3(6) was 
used for studying the K n 2 , and beam background, respectively. See text for more details. 

and photon veto. All the cuts were selected to optimize the background rejection and the 
signal acceptance. Data samples used for this study came from either the 1/3 data (from 
Skim 4 to Skim 6) or the monitor trigger events defined in Section Hi H 31 depending on the 
nature of the selection criteria that were studied. 

1. Single Beam K + Requirements 

The beam cuts were used to identify beam particles scattering either in the beam instru- 
ments or in the target and to ensure a single beam K + particle, making full use of time 
measurements from various sub-detectors, energy loss measurements and pattern recognition 
information in both the B4 and the target as discussed below. 

Beam Times: The tc K , tc n , tsw an d t BM cuts were used to reject extra beam particles 
scattering in the target when one of them agreed with t rs within ±2 ns. Since the incoming 
K + beam intensity was high discriminator dead time was important and the time measure- 
ments used the CCD information in addition to the TDC information from the Ck, CV and 
B4 hodoscope. Fig. [TH shows the single beam K + signal indicated by the Kfi2 monitor 
and the beam background other than the or K n 2 peak events rejected in the kuu(1) 
momentum distribution. 

Energy Loss in BJ^: The tuning of the energy loss cut in the B4 hodoscope used the 
acceptance sample from the K[i2 monitor sample and the rejection sample from the scattered 
7r + 's in the Skim 6 sample. This cut required that the energy loss for an incoming beam 
particle should be consistent with a K + (> 1.1 MeV) as shown in Fig. [T5l 
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FIG. 14: The momentum of charged particles versus tc K — t rs and tc n — t rs . The Kfi2 monitor 
was used to represent the single beam K + events. The data plots indicate beam background 
contamination at beam time (tc K or tc n ) close to t rs in the ttvv{1) trigger sample. The arrows 
indicate the rejected timing regions. The statistics in these plots account for about 0.3% of Nk- 

Delayed Coincidence: The acceptance samples were taken from the peak events in 
the Kfi2 monitor sample, while the rejection sample was from the scattered 7r +, s in the Skim 
6 sample. As shown in Fig. [TBI the distribution of t n — tx for K + decays at rest was an 



exponential, as expected, consistent with the known K + lifetime [18|. The distribution for 
the scattered 7r + 's shows a prompt peak around ns. The delayed coincidence cut required 
tir — tx > 2 ns. When not all sub-detector time measurements were available, the time 
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FIG. 15: The energy loss in the B4 hodoscope for beam K + (solid) and beam ir + (dashed). The 
K +, s were the peak events in the K\x2 monitor sample, while the 7r + 's were the scattered ir +, s 
in the Skim 6 sample. The arrow indicates the cut position. 

resolution was degraded; the delayed coincidence cut was adjusted up to 6 ns delay to take 
into account the resolution. It was expected that this delayed coincidence cut had the same 
rejection power to the background from K +, s decaying in flight. 

Beam Likelihood: The K + stops in the target were required to have energy loss in the 
B4 and target consistent with that expected for the measured K + stopping position. These 
conditions helped to eliminate single beam backgrounds with a scattered n + which did not 
have a consistent path length in the target. The three quantities were combined into a 
likelihood function. Fig. [17] shows the likelihood distributions for signal and background 
using the peak events in the Kfi2 monitor sample and the scattered 7r + 's in the Skim 
6 sample. 

Pileup cut: The pulse shape recorded with the target CCD's provided pileup information. 
The time development of the output signal was fitted with both single- and double-pulse 
assumptions. If the fitted pulse was more likely to be a double pulse and the time of the 
second pulse was coincident with the K + decay time, the event was rejected. 

Pathology Cuts: In the target scintillator, charged particles sometimes underwent nuclear 
interactions or photons had electromagnetic showers. Complicated target patterns could 
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FIG. 16: Distribution of the time difference, t^ — tx measured by the target, for the peak events 
(solid) from the Kfi2 monitor sample and the scattered 7r + 's (dashed) from the Skim 6 sample. 
Beam 7r + 's can be assigned a time tx by the target fiber reconstruction due to mis-identification. 




Beam Likelihood 

FIG. 17: Beam likelihood for K + decays at rest (solid) from the Kfi2 monitor sample and scattered 
7r + 's (dashed) from the Skim 6 sample. The energy loss of K + in the target was required to be 
greater 25 MeV and the number of K + fibers was required to be greater than 2 before evaluating 
the beam likelihood. 
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shift background events into the signal region. Some flaws in the selection criteria were 
identified by inspecting the background events that showed up outside the signal region to 
be described in Section IIII G 81 Below are the situations addressed by pathology cuts. 

• Target 7r + cluster with an identified kink. This was an indication of a hard scattering 
process, which might lead to an incorrect measurement of kinematic quantities. 

• In the B4 the measured energies derived from the ADC and CCD were required to be 
consistent within 1.5 MeV and the measured times derived from the TDC and CCD 
were required to be consistent within 2 ns. Inconsistency in either the times or energies 
was likely due to a second beam particle. 

• Target 7r + fiber next to the decay vertex with energy greater than 3.5 MeV. This 
might indicate a K + fiber being mis-identified as a ir + fiber, causing an incorrect 
measurement of kinematic quantities. 

• Target tt + fiber with energy greater than 5 MeV. This might imply a K n 2 or a radiative 
Kpi event with a photon hiding in the ir + fiber. 

• If a target edge fiber was identified as a K + fiber with a time within 3 ns of the nearby 
IC, the event was rejected. This feature was often an indication of a double beam 
background. 

• Target fibers on the opposite side of the ir + track with respect to the vertex having 
energy greater than 2 MeV within ±3.0 ns of t n . This usually indicated a K n 2 event 
with a photon conversion opposite to the n + . 

• Target tt + track on the opposite side of the 7r + track with respect to the K + decay 
vertex. This usually indicated a double beam or a cosmic ray background in addition 
to the first K + beam particle. 

• Target edge fiber greater than 5 MeV within ±4 ns of t rs . This suggested some photon 
activity in the edge fibers, or a double beam event with the second beam hiding in the 
track. 

• Any activity found in either the UPV or RV within ± 4 ns of t rs . 
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Most pathology cuts were related to target pattern recognition. These backgrounds could 
not be tagged with the usual procedure of background analysis until the outside-the-box 
study given in IIII G 81 

2. Decay 7r + Kinematic Requirements 

Once the ir + kinematic reconstruction routines produced the range-, energy- and 
momentum-related quantities, cuts on the kinematic values (referred to as KIN cuts), were 
identified and divided into several sub-groups which are described below. 

Fiducial Cuts: The n + stopping layer was required to be RS layers 11-18. No charged 
track was allowed to stop in the RSSC layer embedded between the 14 th and the 15 th RS 
layer. Events were also rejected if the stopping layer was 14 with an RSSC hit found in 
the same RS sector or one sector clockwise of the stopping counter. The cosine of the polar 
angle for a charged track was required to be within ±0.5. The z position from the UTC 
track extrapolation to each RS layer was required to be as narrow as ±30 cm to reject 
backgrounds with longer path lengths in RS. The effective UTC fiducial volume was defined 
to be \z\ < 25 cm at the UTC outermost layer. 

Signal Phase Space: The phase space cuts required that the momentum, kinetic energy 
and range of a tt+ track should be in 211 < P < 229 MeV/c, 115 < E < 135 MeV, and 
33 < R < 40 cm. As the resolution of kinematic quantities depended on the azimuthal 
and polar angles of the ir + track, the lower limits of the phase space region were further 
defined by the requirements (referred to as K n2 kinematic cuts) that Pdev = APjap > 2.5, 
Edev = AE/a E > 2.5, and Rdev = AR/a R > 2.75 where the AP, AE and AR were the 
deviation from the K n2 kinematic peak positions, and <jp, a E and op were the corresponding 
resolutions, which were correlated with the azimuthal and polar angles. 

Tracking Quality in Target: Good target tracking relied upon a consistent pattern of hits 
in the target. This required the nearest target K + fiber of the track to the B4 hit position to 
be no more than 2 cm away. The K + decay vertex was required to be nearest to the extreme 
tip of the K + cluster. No more than one fiber gap was allowed between the K + decay vertex 
and the closest ir + fiber. A target track could include either ir + fibers with photons or 
photon fibers mis-identified as the tt + fibers, leading to an incorrect energy measurement 
or event classification if the photon veto also failed. The ir + fibers were therefore examined 
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using likelihoods based on the comparisons of the time, energy, and distance to the track 
between those values expected from the simulated K + —> tx + vv sample and the measured 
values from the Ktt2{\) monitor events. 

Tracking Quality in UTC and RS: At least 4 hits out of 6 cathode foils were required to 
ensure a good measurement of a track in the UTC. Further quality checks were developed 
with respect to the momentum peak, taking account of all the possible circumstances 
that could lead to an incorrect momentum measurement, such as no hit in the outermost 
foil layer, less than 12 layers of anode wire hits, overlapping tracks, or too many wire hits 
in a cluster being excluded from the fit. The cuts were adjusted so that the momentum 
resolution effect did not give a significant contribution to the and fi + background es- 
timates (Section IIII Gp . Another good way to ensure a good UTC measurement was to 
require consistency among UTC, RS and RSSC's. In the x — y plane, this could be achieved 
by requiring consistency between the UTC extrapolation and the positions of the sector 
crossing and the RSSC hits. Similarly, in the <fi — z plane consistency was required between 
the UTC extrapolation and the z positions measured by RSSC and RS. To establish these 
requirements, the /i + 's in Skim 5 were selected with all the other cuts applied, except that 
the maximum momentum cut was not applied and the ir + — > /i + — > e + decay sequence cuts 
(see Section IIII F 3j) were inverted, and the 7r + 's in Skim 6 were selected with all the cuts 
applied, except that the cuts on single beam K + were inverted (Section lHIF II) . These signal 
and background samples were also used in the relevant studies of energy loss in the RS and 
range- momentum consistency in UTC and RS. 

Energy Loss in RS: It was noted that a fi + from the might fake a ir + due to scattering. 
In addition a ir + from the K n2 might fake a signal due to either a photon or an accidental hit 
being hidden in the track counters. The resulting background was removed by comparing 
the energy measurement in each RS layer with the expected value from the range. The 
energy deviation was required to be within ±4<r's for each RS counter as shown in Fig. [181 
In addition to each layer, a probability for energy loss consistency was also calculated with 
all track counters except for the T-counter, the stopping counter and the counters with 
sector crossings. As shown in Fig. [19], the probability cut provided good separation between 
7r + 's and /i + 's. Another way to reject the /i + background was to calculate the likelihood 
from the ratio of the expected (E l expt ) and measured (E l meas ) energy of a track in the z-th 
layer, AE = log E\ xpt — log E l meas , taking into account the Landau tail of the expected energy 
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FIG. 18: Distributions of the maximum energy deviation for 7r + 's (solid) and ^ + 's (dashed) in RS. 
The cut positions are indicated by the arrows. 

distribution. This calculation was done up to and including the layer prior to the stopping 
layer. Fig. |2"U1 shows good ir/ ji separation using the calculated likelihood value. It was noted 
that this cut on the likelihood was correlated with the cuts on the energy deviation and the 
probability given above. 

Range-energy Consistency in IC and Target: Despite the poor resolution of energy mea- 
surements in the IC and target, they could still provide a certain level of particle identifi- 
cation in addition to that from the RS. Signal and background samples used in this study 
were from the Ktt2(1) and Kfi2 monitors. The difference between the Ejc and the expected 
value from Ric was required to be between —5 to 1.75 MeV. The target range {Rt g ) in cm 
and energy (E tg ) in MeV were used to cut events in which R tg > 12 cm, E tg > 28 MeV, 
9.5 x E tg > 28 x R tg or 10 x E tg < 21.5 x (R tg — 2) to reject background with a photon 
hiding along the tt + track in the target. 

Range-momentum Consistency in UTC and RS: This cut was used to check whether the 
range of the charged track was consistent with that for a ir + . The range deviation in RS 
was defined as Xrm = (R rs — Rutc)/o-R, where R utc was the expected range calculated from 
the momentum measured by the UTC with a tt + hypothesis, and an was the uncertainty 
of the measured range as a function of the momentum. The /i +, s and 7r +, s were selected 
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FIG. 19: Distributions of the probability of energy loss consistency in the RS dE/dx measurement 
for 7r + 's (solid) and /i + 's (dashed). The arrow indicates the cut position. 

from Skim 5 and Skim 6 samples, respectively, as used in the study of energy loss in the RS, 
except that the maximum momentum requirement was also applied in Skim 5 to remove 
the K^2 range tail events. The distributions of the range deviation for 7r + 's and /i + 's were 
shown in Fig. [211 Good w/fi separation was observed. 



3. ^ 



e + Decay Sequence 



All cuts for identifying the n + — > /i + — > e + decay sequence were put into a special group 
(referred to as TD cuts). The TD information recorded the pulse shape, providing a tool to 
recognize the decay sequence. As can be seen below, the TD cuts were independent of the 
kinematic reconstruction and photon veto and could be used in the bifurcation study. The 
signature for this decay sequence was: 



Three energy deposits (pulses) corresponding to the 7r + kinetic energy, 7r + 
and /i + — > e + u e u ll decays were found in the stopping counter. 



The kinetic energy of the /i + from ir + — > u^v,,. decay was 4.1 MeV, but due to satura- 
tion the observed energy was about 3 MeV 175]. Since the path length of the /i + was 
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FIG. 20: Likelihood distributions of the RS energy measurement for 7r + 's (solid) and fi +, s (dashed). 
The arrow indicates the cut position. 

~ 1.4 mm in RS, the fraction of /x + exited the stopping counter without depositing 
more than 1 MeV was only ~ 1% of ir + decays. 

• The e + from /i + — > e + i> e i? tl decay has a kinetic energy of E e < 53 MeV. Most of the 
e +, s exited the stopping counter and deposited energy in the other RS counters. 

The three energy deposits from the 7r + — > ji + — > e + decay sequence should be observed by 
the TD's at both ends of the stopping counter. 

For the /i + background, only two pulses due to muon kinetic energy and decay would 
be produced. A fi + could fake a 7r + when an extra pulse was detected in addition to the 
expected two pulses. To suppress the /i + background, two stages of cuts were imposed. 

At the first stage evidence for the 7r + — > decay was sought. The pulse development 
in the stopping counter as recorded by the TD's was fitted with a single- and double-pulse 
hypothesis in an interval of ~4 7r + lifetimes (typically 104 ns). The template shapes used 
in the fit were derived from the average of measured pulses from fi + traversal for each end 
of each RS counter. In addition, a correction was applied to the template shape to take into 
account the change in shape due to propagation along the counter. The parameters of the 
single-pulse fit were the time, the total area of the pulse and a constant corresponding to a 
pedestal of typically 3 TD counts. The parameters of the double-pulse fit were the time of 
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FIG. 21: Distributions of the range deviation in the RS for ir +, s (solid) /i + 's (dashed) tracks. The 
arrow indicates the cut position. 

the first pulse, the time difference of the two pulses, the total pulse area, the fractional area 
of the second pulse and the pedestal. A fit to a triple-pulse hypothesis was attempted if 
evidence for a third pulse was found based on a rudimentary analysis of the TD information 
or if evidence for the tt + — > decay from the double-pulse fit was lacking. The two 

additional parameters in the fit were the time difference of the third pulse with respect to 
the first and the fractional area of the third pulse. The results of the single- and double-pulse 
fit hypotheses for ir + — > y^v^ decay were shown in Fig. [221 Loose requirements were first 
applied on the observed /i + energy 1 < < 14 MeV and on the relative quality of the 
results of the single- and double-pulse fits, R(l) > 1 or R(2) > 2, and R(l) x R(2) > 1 
where the quantities are defined in Table IIVI 

At the second stage, five cuts listed in Table IVl were then applied to suppress the following 
four mechanisms when a fi + — > e + decay could fake the three-pulse ir + — > fi + — > e + decay 
sequence. 

• 7r + time accidental: Accidental activity produced the first pulse, while the /i + — > e + 
decay gave the second and third pulses in the timing sequence. 

• Early /i + decay: The /i + — > e + decay occurred at an early time (< 100 ns), producing 
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FIG. 22: Results of fits with the single- and double-pulse hypotheses to the TD pulse shapes in 
the stopping counter. The solid curve is for the double-pulse hypothesis and the dashed curve is 
for the single-pulse hypothesis. 



Quantity 



Definition or use 



E, 



xl(i) 



Energy of the second pulse 
Time of the second pulse 
X 2 for the n-pulse hypothesis for the counter 
end i {% = 1, 2 =upstream, downstream) 

log 10 (i?(l) • R{2)) Neural net input 
l°gio(Xi(l) • Xi(2)) Neural net input 

dz = z-n — Zn, neural net input component 

dt Time difference between both counter ends 

for the second pulse, neural net input 

TABLE IV: Definitions of quantities determined by the pulse-fitting in the stopping counter. The 
z positions were determined from the energy ratio between the two ends of the stopping counter. 
The z position of the nominal ir and [i pulse was z v or z^ respectively. The term "second" pulse 
identifies the // + -candidate pulse. The neural net is described in the second stage of cuts. 
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in the track counters 






Neural net tt + — ► /i + decay cut 









TABLE V: List of the backgrounds targeted by ir + — > fi + — ► e + decay sequence cuts. 

the second pulse, and accidental activity was identified as the third pulse. 

• /i + time accidental: The fi + — > e + decay made the third pulse, while the second pulse 
was produced by accidental activity occurring between the /i + stop and decay. 

• Tail fluctuation: A fluctuation in the falling edge of the first pulse was identified as 
the second pulse. The decay positron from the fi + decay made the third pulse. 

7r + Time Consistency Cut: This cut suppressed the ir + time accidental background. 
When accidental activity made the first pulse and a charged track made the second pulse 
in the stopping counter, the timing of the first pulse obtained by the TD (t^Ts) wa s not 
coincident with t^s obtained from the other RS counters along the track. Events were 
rejected if \t nt TD — tns\ > 2.5 ns. 

fi + — > e + Decay Requirement: The positron from the /i + — > e + decay (the third pulse) 
generally deposited energy in the stopping counter and other neighboring counters as de- 
picted in Fig. [23j The positron finding started by looking for a cluster of TDC hits in the RS 
counters in the region within ±1 sector and ±2 layers of the stopping counter. The cluster 
should include the stopping counter and an additional hit in the same sector as the stopping 
counter. Candidates for the positron were found when the average time of the hits in the 
cluster was within ±2.4 ns of the TDC time (the third pulse) in the stopping counter. The 
z positions of the hits in the cluster, obtained from the end-to-end time differences of the 
hits, were also required to be consistent with the z position in the stopping counter. If the 
candidate was due to a track that passed through the stopping counter, then hits might be 
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FIG. 23: Schematic view of signal (left) and background (right) of the /i + — > e + decay. Each 
rectangle represents an RS counter. The central counter represents the stopping counter and the 
shaded rectangles represent hit counters. The arrow indicates the possible direction of the positron 
or charged track producing the hits. 



found on both sides of the stopping counter (Fig. 1231) . The early /i + decay background was 
removed by requiring that the second pulse from the TD pulse fitting agreed with the time 
of the cluster. 

Cut on Accidental Activity: Accidental activity in the stopping counter was frequently 
associated with activity in other RS counters as well as the BV, BVL and EC. Hence, 
events with activity coincident with the second pulse in the stopping counter were targeted 
for rejection. The time windows and energy thresholds for the various subsystems in RS, 
BV, BVL and EC were optimized in order to have the highest rejection power at a given 
acceptance value of 94% for this cut. Events were rejected if the energy sum of the hits 
within a time window in any of the subsystems was greater than the threshold. There was 
also a kind of accidental activity that overlapped the charged track and made a second 
pulse in the stopping counter. To reject this accidental background, fits were performed 
to a double-pulse hypothesis in the two RS counters along the track prior to the stopping 
counter. If the time of the fitted second pulse was within ±5 ns of the second pulse in the 
stopping counter, the fitted energy of the second pulse was greater than 1 MeV and the \ 2 
ratios of the single- to the double-pulse fit hypotheses were greater than 4, then the event 
was rejected. 

Neural Net 7r + — > fi + Decay Cut: The tail fluctuation background mimicking the energy 
deposit for a fi + at the falling edge of the 7r + -induced pulse was characterized by a small 
decay time and a low pulse area in the second pulse. The variables shown in Fig. [24] and 
described in Table HVl differed for events induced by tt + — > and /i + — > e + v e P (U decays. 
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FIG. 24: Distributions of the input variables for the Neural Net function in n + — > [i^v^ induced 
(solid) and (i + — > e + v e D^ induced (dashed) decays. The energy of the second pulse (a), time of 
the second pulse (b), log of the product of the single-pulse fit x 2 's for both ends (c), z position 
difference between the first and second pulses obtained from the energy ratio of both ends (d), log 
of the product of the x 2 ratios of single- to double-pulse fits for both ends (e) and time difference 
of the second pulses in both ends (f). 

Application of a fixed cut to each variable would cause a non-negligible acceptance loss. In 
order to achieve a higher acceptance at the same rejection as the fixed cuts, a Neural Network 
(NN) technique was adopted. The NN function was derived via a Multi-La yer Perception 
program incorporated in the library of Physics Analysis Workstation (PAW) 76[ . To create 
the NN function, the scattered 7r + 's and range tail events in the 7rz/i/(l) trigger, which 
passed all other 7r + — > fi + — > e + decay sequence cuts, were used as signal and background 
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FIG. 25: Distributions of the outputs of the NN function for 7r + (solid) and [i + (dashed) events. 
Events with an output of the NN function less than 0.76 were rejected. 

samples, respectively. A 5-variable NN function was obtained using the six variables shown 
in Fig. [23] with the differences in z position and time combined to create a single input 
variable, 

X 2 (z,t) = (dz/a dz ) 2 + (dt/a dt ) 2 , (24) 

where the z position and the time were obtained from the energy ratio of both ends, and 
time difference of the second pulses in both ends. The resolution of dz (dt) was a dz 
Fig. [25] shows the distributions of the output of the NN function for 7r + 's and /i + 's. The 
rejection of the NN ir + — > /i + decay cut as a function of the acceptance is shown in Fig. [26] 
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FIG. 26: Rejection of the NN ir + — > /x + decay cut as a function of the acceptance. The intersection 
of the vertical and horizontal dashed lines shows the rejection and acceptance at the nominal cut 
position. 



4- Photon Veto 

To achieve the background level much less than one event, the total 7r° rejection was 
required to be of order of 10 6 . A rejection factor of ~ 10 4 was already achieved online, leaving 
a further ~ 10 2 rejection factor to be achieved by the offline analysis. The corresponding 
photon veto cuts (referred to as PV cuts) were used to identify the photon activities detected 
by all the PV counters. 

A search for the photons coincident with the track time was performed in the subsystems 
of the BV, BVL, RS, EC, target, IC, VC, CO and /iCO. The timing resolution of each photon 
detection system was a key ingredient in determining the time window for the PV due to the 
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false veto rate in the high rate experimental environment. Fig. [271 shows the measurements 
of the timing resolution for each system as a function of the visible energy. Since the EC 
consisted of 4 rings and the inner ring (ringl) was exposed to high accidental rates from 
the beam, the corresponding time resolution was worse. The time window and energy 
thresholds in each subsystem were optimized by adjusting the cut positions to maximize 
rejection for a given acceptance. The rejection sample was from Skim 4 with K^i events, 
while the acceptance sample was from the K[i2 monitor. To ensure the events were from the 
K n2 peak, the measured range, energy and momentum of the charged track were required 
to be within three standard deviations of the nominal values. For the K^, the momentum 
of the charged track was required to be within three standard deviations of the nominal 
value, and the range was required to be longer than 37 cm to remove any event which could 
contain photon(s). The selected K^i and events were also required to pass the KIN 
cuts and the beam cuts. 

The optimization process [3] started from the initial set of cut parameters. For a new 
set of parameters, the rejection and acceptance were re-measured. Only one subsystem's 
cut parameters were varied at a time. If the rejection increased without losing acceptance 
or the acceptance increased without losing rejection, the set of parameters was regarded as 
a good input for the iteration. More preferable cases occurred when both the rejection and 
acceptance were improved. Fig. [281 illustrates this optimization process for the photon veto. 
The optimization process continued until no more gain was obtained in rejection without 
losing acceptance. The boundary point was measured at every given acceptance position. 
As a result, a profile curve, which gave the maximum achievable rejection, was obtained. 

Fig. [29] shows the offline rejection of the photon veto cuts against the K n2 background 
as a function of acceptance. The time window and energy threshold for each category are 
tabulated in Table IVI1 As a reminder, the total offline photon veto is not the simple product 
of the rejections listed in Table [VI] due to mutual correlations. 

G. Background Evaluation 

To have an unbiased result, the signal region was always masked until all the background 
evaluation studies were completed. As described in Section IIII Al the stopped K + decay 
background and the beam background were subdivided into the following categories: 
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FIG. 27: The timing resolution as a function of visible energy in various PV counters. 



• K n2 background, 

• K^ 2 range tail background, 

• fi + band backgrounds, 

• single beam background, 

• double beam K + — K + background, 

• double beam K + — n + background, and 
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Category Time window (ns) Energy threshold (MeV) Rejection 



BV 


±4.50 


0.20 


14.556 


BVL 


±2.00 


0.00 


1.247 


RS 


±1.50 


3.80 


3.329 


EC 


±2.25 


3.80 


2.296 


Inner EC 


±1.75 


1.00 


1.186 


Target 


±1.30 


4.79 


1.493 


IC 


±2.25 


0.40 


1.321 


vc 


±2.00 


1.60 


1.054 


CO 


±1.25 


0.80 


1.052 


//CO 


±1.50 


2.80 


1.003 



TABLE VI: Time window and offline energy threshold for each category of the photon veto cuts. 
Also listed are individual rejection values contributed by each sub-system. 
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FIG. 28: Illustration of the optimization process to determine the photon veto parameters. 
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FIG. 29: Offline rejection of the photon veto cuts against the K n 2 background as a function of 
the acceptance. The crossing point of the vertical and horizontal lines shows the rejection and 
acceptance at the cut position. 

• CEX background. 

Except for CEX, all background levels were estimated using the data by means of the 
bifurcation method. Table IVHI gives the bifurcation cuts, the data stream categories and 
the results for the 1/3 data sample. The background levels given in Table IVTT1 must be 
scaled by a factor of 3 to obtain estimates for the full sample. The 1/3 sample was studied 
first in order to tune and optimize the cuts. Then the 2/3 sample was used to give the final 
background estimates. When the normalization branch contained only one or few events 
in region B, a second bifurcation analysis was performed in this branch to improve the 
statistics. Details of these procedures are given in the following sections. 
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CUT1 


CUT2 Category 


B (C + D)/C 


BC/D 


K o 


PV 


KIN 


S\lriTn ZL 


0.39 ±0.11 85.2 ±2.5 


0.0046 ± 0.0013 




TD 


KIN 


Skim 5 


1.81 ±0.16 (4± 1) x 10 2 


0.0041 ±0.0011 


band 


TD 


KIN 


Skim 5 


2.37 ±0.74 (4± 1) x 10 2 


0.0053 ± 0.0021 


Single-beam 


DC 


B4 


Skim 6 


8 (7 ± 4) x 10 3 0.0011 ± 0.0007 


Beam K + — 


K+ BWPC B4 


Skim 6 


0.04 ±0.04 117 ±37 


0.0003 ± 0.0003 


Beam K + — 


vr+ BWPC B4 


Skim 6 


0.26 ±0.11 (7 ±4) x 10 3 


< 0.0001 



TABLE VII: Results of the bifurcation analyses for the backgrounds in the 1/3 sample only. All the 
results for the region B were from the second bifurcations except for the single-beam background as 
explained in the text. Details for CUT1 and CUT2 are described in the text. Errors are statistical 
only. 

1. K n 2 Background 

A K^2 decay event should have a charged track with a monochromatic momentum, range 
and energy plus two photons. Experimentally, if a K n2 event appeared in the signal region, 
the photons from the 7r° decay must have escaped detection and the K n2 kinematics must 
have been distorted by scattering or resolution effects as well. 

In the background study, the two bifurcation cuts were chosen as the PV cuts (CUT1) and 
the signal phase space cuts in the KIN cuts (CUT2), since both of these could independently 
give powerful rejection of the K^i background. In order to remove the contamination from 
/i + and beam events, the bifurcation analysis sample was selected from the Skim 4 sample 
by applying the TD cuts, the beam cuts and the KIN cuts other than the phase space cuts. 
In the rejection branch, K n2 events were selected by inverting the signal phase space cut 
in the KIN cuts (CUT2), giving 95,797 events for the region C+D. The PV cuts (CUT1) 
were then applied to the remaining K n 2 events, leaving 1,124 events for the region C. In the 
normalization branch, the K n 2 events with photon activity were selected by inverting the 
PV cut (CUT1). The signal phase space cuts in the KIN cuts (CUT2) were applied to the 
above selected K n2 sample, resulting in no events (B = 0) left in the normalization branch. 

To deal with the above situation and give a non-zero events in the normalization branch, 
another (second) bifurcation analysis was performed by separating the kinematic cuts 
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into Edev (CUT1) and Rdev+Pdev cuts (CUT2), since the E measurement was almost 
independent of R and P measurements. These two bifurcation cuts were applied sequentially 
to the selected K n2 sample in the normalization branch. In this second bifurcation study, 
the lower boundary cuts on E, R and R were removed. Changing the cut positions on the 
Rdev, Pdev and Edev gave the number of events in the normalization branch (Nrp) and the 
rejection branch (Re) in the second bifurcation analysis, which were then used to calculate 
the expected number of events in the normalization branch by means of B = Nrp/(Re — 1). 
This second bifurcation method gave smaller uncertainty, providing a way to optimize the 
phase space cut positions to reject the K v2 background with less acceptance loss. The 
expected number of events from the second bifurcation analysis was found to be ~ 10% 
less than the observed number from the first bifurcation analysis in the region B+D. This 
was due to the small level of correlation between R and E when estimating the ir + range 
in the stopping layer from the measured energy in the stopping layer. The cut positions 
were chosen at Rdev > 2.75, Edev > 2.5 and Pdev > 2.5, in order to reduce the expected 
K n 2 background level to about 0.01 events level as shown in Table |VlTl Fig. 1301 shows the 
expected K n 2 background from the second bifurcation in the 1/3 normalization branch as 
a function of the Edev cut position. The acceptance was measured using the simulated 
K + — * ir + uu sample. 

The expected number of events from the second bifurcation was used to give the result 
for the normalization branch, which was then used to give the estimated background level 
for the 1/3 sample as shown in Table IVHl 

2. fi + Background 

The n + background consisted of the range tail events and the /i + band events as 
indicated in Fig. El These /i + 's lost energy and eventually came to rest in RS and could 
migrate into the signal region through resolution effects if the TD cuts failed. 

The two bifurcation cuts with the most powerful rejection were the TD cuts (CUT1) 
and the selected KIN cuts (CUT2), which excluded the fiducial cuts, the lower boundary in 
the signal phase space cuts, the cuts on the tracking quality in target and the cuts on the 
range-energy consistency in IC and target as detailed in Section IIIIF2[ The fi + background 
sample used for the bifurcation study was selected from the Skim 5 sample by applying the 
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rm 



FIG. 30: The expected kinematic background events in the 1/3 normalization branch (histogram 
with left axis) as a function of K n 2 kinematic cut (Edev) position for the K n 2 background (a), 
maximum momentum cut position for the Kai range tail background (b) and range-momentum 
cut (xrm) position for the fi + band background (c). Also shown is the relative acceptance change 
as the cuts (curve with right axis). 



PV cuts, the beam cuts and the KIN cuts not used in this bifurcation study, in order to 
remove the K n 2 and beam backgrounds. In the rejection branch, there were 7,119 events 



in the region C+D when inverting the selected KIN cuts (CUT2). The TD cuts (CUT1) were 
applied to the /i + events, leaving 16 events in the region C. In the normalization branch, 



the n + events were obtained by inverting the TD cuts (CUT1). The selected KIN cuts 
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(CUT2) were applied to the above selected /x + sample, leaving only one event (B = 1) in 
the normalization branch. 

The result from the bifurcation analysis given above had a large statistical uncertainty, 
since only one events remained in the normalization branch. The K^ 2 range tail and the /x + 
band events were the only two possible fi + backgrounds in the signal region. The origins 
of these two backgrounds were due to the momentum and range resolution effects. Once 
these resolution effects were known in the signal region, a better estimate of corresponding 
events B in the normalization branch could be obtained. In this analysis, the events with 
P > 225 MeV/c were regarded as the K^ 2 range tail background. 

The momentum resolution effects could be well described by using the momentum dis- 
tribution from the K^ 2 peak events. In order to enhance the number of events in the nor- 
malization branch, the RS energy loss cuts were removed from the KIN cuts. Also removed 
were the upper boundary cuts on the R, E and P. The K^ 2 peak events were selected by 
requiring R > 50 cm. The selected K^ 2 peak events in the 7rz^(l) trigger were found to have 
longer range in the target, leading to a bias of 0.5 MeV/c higher momentum measurement 
when applying a ir + hypothesis to the contribution from the energy loss in the target. After 
subtracting this bias for each K^ 2 peak event, the momentum distribution was seen to be 
in good agreement with that from the K^ 2 range tail events and with more statistics in 
the signal region defined below 229 MeV/c. Normalizing the number of peak events to 
that of the K^ 2 range tail events observed in the region B+D gave the expected number of 
events B in the normalization branch. The expected range tail background is given in 
Table IVHI Fig. [30] shows the number of K^ 2 range tail background events as a function of 
the maximum momentum cut position in the normalization branch. 

The /i + band background came from the RP resolution effects, which were well described 
by the range deviation (Xrm) in RS. Enhancing the number of events in the study of the 
RP resolution effects was achieved by removing the PV cuts and the RS energy loss cuts 
in the normalization branch sample, since they were not correlated with the range and 
momentum measurements. Within the statistical uncertainty, both distributions were seen 
to be consistent except that the distribution without the above requirements gave higher 
statistics in the normalization branch (Fig. 130]) . The estimated background level for the /i + 
band is given in Table IVHI 
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3. Single Beam Background 



In the study of single beam background, the two bifurcation cuts chosen were the offline 
delayed coincidence (DC) cuts (CUT1) and the B4 energy loss cut (CUT2). The offline DC 
cuts were from the precise offline time measurements from the beam instrumentation, the 
target, the IC and the RS (They were not the ones used in the trigger). The events were 
selected from the Skim 6 sample by applying the PV cuts, the KIN cuts, the TD cuts, and 
the beam cuts except for the DC cuts and the B4 energy loss cut. In the normalization 
branch, the single beam events were selected by inverting the DC cuts (CUT1). The the 
B4 energy loss cut (CUT2) was applied to the above selected single beam sample, leaving 8 
events in the region B for the normalization branch. In the rejection branch, 29,100 single 
beam events in the region C+D were selected by inverting the B4 energy loss cut (CUT2). 
The DC cuts were applied to these selected events, resulting in 4 events in the region C. The 
rejection factor was applied for both the ir + scattering events and the K + decay-in-flight 
events, since there was no reason to have different rejections for these. The background 
estimate for the single beam background is given in Table IVHI 

4- Double Beam Background 

As already defined in Section IIII Al the double beam background could be due to a 
K + — K + event or a K + — 7r + event. For this double beam background, the DC cuts were 
insufficient to remove the double beam background, but the time difference between the 
beam instrumentation (Ck, C n , BWPC's and B4) and the 7r + track was a good indicator. 
Another independent way was to use the target pattern to identify extra particles other than 
the initial K + hit. 

The two bifurcation cuts chosen were the BWPC timing cuts (CUT1) and the B4 timing 
cuts (CUT2). The events were selected from the Skim 6 sample by applying the PV cuts, 
the KIN cuts, the TD cuts and the beam cuts except for the BWPC and the B4 timing 
cuts. These cuts removed the K^, ^ + and single beam backgrounds. In the normalization 
branch, there were no events left if the B4 timing cuts (CUT2) were applied. In order to give 
a more precise estimate, the second bifurcation was adopted. Time measurements including 
the trailing edge TDC time from the Ck and C n were used to select K + — K + events and 
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K + — 7r + , separately. The bifurcation analyses were then performed using the B4 timing 
cuts and the target pattern recognition cuts. Results are given in Table IVHl In the rejection 
branch, the double beam background events were selected by inverting the B4 timing cuts 
(CUT2). The K + — K + and K + — ir + events were tagged by the Ck and C n , separately; 
this resulted in 1,170 and 22,150 events for both cases. Applying the BWPC timing cuts 
(CUT1) resulted in 10 and 3 events observed for the K + — K + and K + — tt + backgrounds, 
respectively. The resulting rejections and background estimates are given in Table IVTT1 The 
K + — K + background was found to dominate the double beam background. 



5. Charge Exchange Background 

Since there was no reliable way to isolate the CEX events from the itvv{l) trigger data, 
the background study could only rely on the Monte Carlo simulation. The CEX simulation 
needed a number of inputs, such as the CEX re-generation rate as a function of K\ energy, 
the K\ decay vertex and the K\ momentum, all of which could only be obtained from 
the real data. A special CEX monitor trigger as described in Section III H 31 was used for 
collecting data with two charged tracks from the K s decay. The A^'s were reconstructed in 
the 7r + 7r~ decay mode and used to measure the K s production rate, momentum spectrum, 
the B4 hit information, decay vertex distribution, and pattern of target A + fibers' time and 
energy. Since a K° decays approximately equally to K\ and K s states, the measured rate 
of K s decays can be used to obtain the K\ production rate in the target 

N K ° S 

Rk1 S e K o-A P yB(Kl Jn+n-)-N K /PS' (25) 
= 2.73 x 10" 5 , (26) 

where the quantities used in this calculation are summarized in Table IVIII1 

Both K*l — > TT + fi~u /1 (Kt^) and K\ — > n + e~u e (K® s ) decays were simulated to estimate 
the CEX background level. Each decay mode was generated with the amount of K° L decays 
equivalent to 7.604 x 10 15 A'£'s. These Monte Carlo events were passed through all of the 
selection criteria, except for the 7r + — > fi + — ► e + decay sequence cuts and the beam cuts that 
were not related to target quantities. There were 56 K® 3 events and 21 A° 3 events surviving 
in the signal region. Therefore, the expected CEX events can be estimated by means of 

7§ 



Ncex = (jV£fl + x ^ X Fa, 
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Description 



Parameter Values 



Number of selected Kg events N K o 8,086 

Kg selection efficiency e^o 0.138 

Photon Veto efficiency Apy 0.680 

K° s -> vr+vr^ Branching Ratio B{K° S -» ir + ir~) 0.686 

Number of K + Triggers (10 12 ) N K 1.77 

Prescaling Factor PS 384 

K° L Production Rate R K o 2.73 x 10" 5 

TABLE VIII: K^ production rate and the quantities that are used to estimate it. 

1 77 x 10 12 

= (56 + 21) x — — x 0.250 

v ; 7.604 x 10 15 

= 0.00448 ±0.00051(stat), (27) 

where Np a f s and Npjfs are the numbers of K® 3 and K® 3 events surviving all the cuts, N^ c 
was the total exposure of generated K +, s and F acc was the acceptance of the TD and beam 
cuts that were not applied to in the simulation. 



6. Initial Background Evaluated from 1/3 Sample 

The initial total background evaluation based on the 1/3 sample was 0.05 ± O.OIsm 
events, which came from the results in Table IVlTI and Equation (j2"Tj) . Given this relatively 
low background level, the signal region was expanded to gain acceptance at the cost of 
more background. In addition to the total background level, the bifurcation analyses gave 
the corresponding predicted background functions for the TD, the PV and the kinematic 
cuts, with which the expected background level and relative acceptance change for a given 
cut position were obtained. These functions were used in optimizing the selection criteria, 
studying the correlation of two bifurcation cuts and determining the branching ratio as well. 



7. Optimization of Signal Region 

The distributions of signal and background in the cut space were described well by the 
predicted background functions and were used to enlarge the signal region. The cuts to 
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be loosened were the NN tt + — > fi + — > e + decay cut in the TD cuts (Fig. [261) . the PV cuts 
(Fig. [291) and the K^i kinematic cuts (Fig. 130]) . Loosening the cuts for the beam backgrounds 
and the other TD and KIN cuts did not provide much acceptance gain. Simultaneously 
loosening the cut positions of the NN n + — > /i + — > e + decay, PV and K n2 kinematic 
cuts could increase the background levels to an unacceptable level. Instead, when one of 
the three cuts was loosened, the cut positions of the other two were kept unchanged. The 
revised (extended) signal region consisted of the standard signal region plus three extensions. 
Hereafter, the revised and standard signal regions were referred to as the extended signal 
region and the standard region, respectively. The regions created by loosening the NN 
7r + — > /i + — > e + decay cut, the PV cuts, and the K n2 kinematic cuts were referred to as 
7r + — > fi + — ► e + extended, PV extended, and K n2 kinematic extended regions, respectively. 
All of which included the standard regions. The loosening factors for these three cuts, 
fwfie = 4.2, fpv = 4.0 and fx^ 2 = 9-5, meant the corresponding background increased by 
the same factor. The total estimated acceptance gain by enlarging the signal region was 
31% (12% from the NN tt+ -> /i+ -> e+ decay cut, 7% from the PV cuts, and 12% from the 
K n2 kinematic cuts). 

8. Correlation and Single Cut Failure Study 

The bifurcation procedure described above assumed that the two bifurcation cuts were 
not correlated. This assumption was tested by comparing the predicted and observed rates 
near but outside the signal region. 

A schematic representation of the region near, but outside, the signal region is shown in 
Fig. EH A near region outside the signal region (A' without the black region in Fig. |3T|) 
was defined by loosening two bifurcation cuts (CUT1 and CUT2) simultaneously by the 
fixed factors a and b, respectively. If the predicted background functions associated with 
the CUT1 and CUT2 were correct, the expected number of background events in the near 
region (BG ! ) was estimated by the bifurcation method as 

BG' = B'C'/D' - BC/D (28) 

If a deviation was seen between the observed and the predicted numbers in the near region, 
then a correlation between the bifurcation cuts could be indicated and the background 
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FIG. 31: Pictorial explanation of the correlation study using the events near, but outside the signal 
region. The vertical axis "a" and the horizontal axis "6" are for the loosening factors, while the 
lxl black region represents the standard signal region. 



estimate might be unreliable. The observation was performed in the same way as the 
bifurcation method used in the background estimation in the signal region, except that two 
bifurcation cuts were loosened by factors of a and b. 

The results of the correlation study for the K n2 , range tail and /i + band backgrounds 
are summarized in Table 1X1 for the 2/3 data sample. Good agreement was found between 
the observed and the predicted number of events. The test results were obtained from the 
comparisons between the observed and the predicted numbers of events using the predicted 
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background function method. 

In addition to the study of correlation between the cuts used for the bifurcation method, 
events that passed all except for a single cut were examined to determine if each cut operated 
as designed for the appropriate background mechanism. Such a study provided a way to 
discover any new type of background or potential analysis flaw. In the 1/3 sample, six of the 
eight events that failed a single cut only were far from the cut position, while the other two 
events showed potential analysis flaws. The first flaw would artificially increase the measured 
range and momentum of 7r + 's from K v2 decays that exited the upstream end of the target 
through the gap between the front face of the target and the B4 hodoscope. Additional cuts 
with minimal acceptance loss were devised to eliminate such events. The second revealed a 
possible correlation between the PV cuts and the KIN cuts when using the tt + polar angle 
(6) as a reference to exclude the accidental hits in the opposite site of the PV counters. 
The corresponding calculation used in the PV cuts was subsequently removed. All the cuts 
designed at this stage were referred to as pathology cuts as described in Section IIIIF II 

9. Final Background Evaluated from 2/3 Sample 

Evaluation of the final background levels came from the 2/3 sample. To get the values 
for the extended box, the corresponding values were scaled by the loosening factors given in 
Section fill G 71 It was noted that loosening the NN 7r + — > fi + — > e + decay cut, the PV cuts, 
and the K v2 kinematic cuts could lead to a small change in the beam background levels. 
Therefore, the beam backgrounds were re-estimated in the extended signal region. The TD 
rejection and PV rejection were measured to be 445 ± 111 and 84.3 ± 1.2 with the 2/3 sample 
in good agreement with those obtained from the 1/3 sample given in Table IVHl The final 
background estimates are summarized in Table IIXL 

The total background level in the extended signal region was estimated to be 0.30±0.03 s t a t 
events, which was dominated by the K n2 background contribution. As the background distri- 
bution was not uniform in the signal region, the predicted background functions obtained in 
the background study were exploited to interpret any possible candidate events observed and 
to give a proper branching ratio measurement by using the likelihood technique described 
in Section HVl 
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Background Standard Extended 

K n2 0.019±0.004 0.216±0.023 

K,j2 range tail 0.010±0.001 0.044±0.005 

band 0.005±0.002 0.024±0.010 

Single-beam 0.004±0.002 0.006±0.002 

Double-beam 0.003±0.002 0.003±0.002 

CEX 0.004±0.001 0.005±0.001 

Total 0.05±0.01 0.30±0.03 

TABLE IX: Total background level and the contribution from each background source for both 
standard and extended signal regions as estimated from the 2/3 data samples. The errors are 
statistical. 

10. Systematic Uncertainty 

Systematic uncertainty in the background estimates arose from the possible correlation 
between the two bifurcation cuts. The correlation was investigated using the 2/3 sample 
and the results are given in Table |2D In addition, the ratios of observations over predictions 
were used to quantify the degree of consistency and were found to be consistent with unity 
within a relative uncertainty of 15%, confirming that the background estimations obtained 
with the bifurcation method were reliable. 



H. Acceptance and Sensitivity 

To reduce the estimated background level to less than one event in the signal region, this 
analysis utilized many selection criteria. The corresponding acceptances for these selection 
criteria were estimated directly from the data when possible, by splitting them into compo- 
nents that could be measured separately using the monitor trigger data or the Monte Carlo 
simulation. The latter gave the estimates on the decay phase space, the trigger efficiency 
and the nuclear interaction effects. 
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0.3±0.1 1.3±0.4 


3.2±0.9 


5.2±1.5 


10.4±2.8 


Observation 


1 


1 


4 


5 
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TABLE X: Results of a correlation study for the K n 2 (top), range tail (middle) and u + band 
(bottom) backgrounds in the 2/3 sample. The errors in the predictions are statistical uncertainties. 

1. Acceptance Factors from K^ 2 Events 

Since the K^ 2 events have the same features as the signal regarding the K + beam, the 
charged track and the event topology, the acceptances associated with the relevant cuts as 
listed in Table IXII were directly measured using the Kfj2 monitor trigger data. Below are 
the details of these measurements. 

Tracking in RS: To measure the acceptance of RS tracking, additional requirements ( "set- 
up cuts" ) were applied to K/i2 monitor data to ensure a good track in the RS without using 
RS measurements. The events were required to have an IC hit with tic ~ tc K > 5 ns, at 
least 1.2 MeV energy loss in B4, and successful tracking in both the UTC and target. All 
surviving K^ 2 events were examined for consistency with RS tracking. The acceptance of 
the RS tracking cuts is given in the second row of Table IXII 

Tracking in UTC and Target: Because there should be no photon activity for peak 
events, the KyJl monitor was an ideal sample to measure the acceptance of target pattern 
recognition criteria. The sample was taken from events surviving the RS tracking cuts 
discussed above. To eliminate possible beam background contamination, this sample was 
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Cut Acceptance 

Tracking in RS 0.99996±0. 00001 
Tracking in UTC and target 0.99568±0. 00010 

Beam selection criteria 0.50779±0. 00074 

Photon veto 0.76784±0.00218 

Track stop in RSSC 0.98195±0.00015 

Muon veto in RS 0.99591±0.00057 
A K ^ 0.3796±0.0013 

TABLE XI: Acceptances of the K + — ► ix^vv selection cuts measured from the K^ 2 monitor trigger 
data. The acceptance of beam cuts does not include those using the energy measurement in the 
target. The errors are statistical. 

also required to meet the timing requirements on the beam instruments. A 5 ns timing 
consistency was also required between tic and t rs . A subset of the PV cuts was applied to 
suppress possible K^ 2l contamination. The BV and BVL elements of the PV cuts were not 
applied to avoid self-vetoing by long range K^ 2 events. Events surviving the set-up cuts were 
then checked with the UTC and target requirements except for those involved ir + energy 
and range measurement in the target, giving a measurement of the acceptance of tracking 
in UTC and target (Table IXII) . 

Beam selection criteria: The K^ 2 events were chosen from those satisfying the require- 
ments on tracking in the UTC and target described above but without the timing require- 
ments on the beam instruments. To suppress beam background contamination, the mo- 
mentum deviation was required to be within two standard deviations of the K^ 2 peak with 
|cos#| < 0.5. Also required was that there should be no discernible scattering of tracks in the 
RS. The remaining events then passed through the beam cuts except for the pathology cuts 
using 7r + energy and range measurement information described in Section IIHF 11 providing 
a measurement of the corresponding beam selection acceptance. Because K^ 2 events were 
simple single tracks, the efficiency of the DC trigger was also measured and included in the 
acceptance of beam selection criteria (Table IXll) . 

Photon veto: The acceptance of the PV cuts included contributions from both the online 
and offline PV. Ideally, K^ 2 events should not contain any photons, and could thus be used 
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to measure the acceptance loss due to the PV cuts. In the first step of this procedure, 
the selection criteria were applied to remove possible beam backgrounds. However, it was 
noted that some /i + 's could penetrate the whole RS and reach the PV counters, resulting in 
a time-coincident PV hit and therefore an over-counting in the acceptance loss due to the 
PV cuts. To avoid this problem, the selected sample was further required to have the 
stopping layers prior to the 19 th RS layer. The PV cuts were then applied to the surviving 
K^2 events, yielding a measurement of the acceptance loss due to the application of the PV 
cuts. Since the K/jl2(1) trigger did not have an online PV cut applied, this acceptance factor 
also included the contribution from the online PV (Table IXllh 

Track stop in RSSC: This cut was classified into the fiducial cuts in Section IIII F 21 and 
aimed at vetoing possible associated photon activity detected by the RSSC, even though it 
was not included in the PV cuts. Using the Kfi2 monitor events and applying the above 
cut to those with hits in the second layer of RSSC gave a measurement of the corresponding 
acceptance value (Table IXlT) . 

Muon veto in RS: The ~kvv(1) trigger condition 19 c t was also called as a muon veto 
in RS. This trigger requirement could result in acceptance losses when an accidental hit 
happened in the 19 th layer along with an otherwise good signal candidate event. This loss 
was measured with the K[i2 monitor events, which were selected by requiring the stopping 
layer to be RS layer 17 and the range to be longer than 40 cm, in addition to applying all 
the cuts used in the above studies except for the momentum cut. In this selected sample, 
the online trigger condition 19^ was checked, giving a measurement of the acceptance for 
the muon veto in RS (Table IXll) . 

2. Acceptance Factors from K n 2 Events 

The Kn2(l) monitor trigger data were used to measure the acceptances for the pathology 
beam cuts involving the ir + energy and range measurement and the KIN cuts involving the 
range-energy consistency in IC and target. This was a complement to the measurement of 
acceptance factors from the events. The K n2 events were selected with all cuts applied 
except for those to be measured. To ensure good K n 2 events, the momentum, range and 
energy were required to be within two standard deviations of the K n2 peak positions and 
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observation of a 7r° — > 77 decay was required. The result was 

A Kn2 = 0.8785 ± 0.0029 stat . (29) 

3. Kinematic Acceptance from Beam ir + Events 

The TT SC at monitor trigger data provided a pure 7r + sample to measure the KIN cuts 
related to the particle type: the cut on 7r + stopping layer in the fiducial cuts, the cuts on 
tracking quality in UTC and RS, the cuts energy loss in RS and the cut on range- momentum 
consistency in UTC and RS. The events were required to pass the Pass 1 cuts and the TD 
cuts. The K + selection criteria in the beam instruments were inverted to select beam 7r + 's. 
The tic was required to be within ±5 ns of t rs . The signal phase space cuts were additionally 
applied to select the events. The acceptance was measured to be 

A nscat = 0.6161 ± 0.0085 stat ± 0.0189^. (30) 

Since these 7r +, s came from the beam 7r +, s scattering in the target and not from the K + 
decays at rest, classification of the K + fibers and tt + fibers could be complicated because of 
their nearly coincident times and differences in fiber energy deposits of scattered 7r +, s and K + 
decays at rest. Both of these features would result in more uncertainties in the momentum, 
energy and range measurements in these 7r + events. Systematic uncertainties were therefore 
investigated by loosening or tightening the signal phase space cuts by ±1 standard deviation. 
The corresponding variation in acceptance was treated as the systematic uncertainty. 

4- vr + — ► n + — ► e + Decay Acceptance from Beam ir + Events 

The acceptance of 7r + — > /i + — > e + decay sequence cuts was measured by using the 7i scat 
monitor trigger data. This acceptance measurement included the online and offline (TD) 
7r + identification cuts. The online ones included LI A and LI. 2 in trigger. The sample was 
selected using the same cuts as those used in measuring A nBcat except for the cuts to be 
measured here. It should be noted that this measurement included the acceptance loss due 
to the 7r + absorption and tt + decay- in-flight. This loss was estimated to be 1.4% using Monte 
Carlo in Section IIII H 51 and should be corrected to remove the double counting problem in 
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the acceptance. The final acceptance of ir + — > /i + — > e + decay sequence cuts was given below 

A n ^ e = 0.3523 ± 0.0077 stflt ± 0.0067^. (31) 

Since the ir + — > fi + — > e + decay sequence cuts might be correlated to particle identification 
KIN cuts when using information from the RS, the effect on acceptance was investigated 
with and without the RS energy loss cuts. The observed about 2% variation on the relative 
acceptance was assigned as the systematic uncertainty. 

5. Acceptance Factors from Monte Carlo Simulation 

Monte Carlo simulations of the K + — > tt + vv were used to evaluate the trigger acceptance, 
the phase space acceptances and the acceptance loss due to 7r + absorption, decay- in-flight, 
and nuclear interaction, which could not be measured directly by the monitor trigger data. 
The K + — > tx + vv Monte Carlo samples were generated with and without including the nu- 
clear interaction. 

Trigger Requirements: All the trigger requirements as described in Section III HI were 
simulated by Monte Carlo except for the DC, Ll.l and LI. 2, which were already measured 
using the K^ 2 and 7r scat monitor trigger data. From the Monte Carlo without including the 
nuclear interaction, the acceptance for the trigger requirements was 

A trig = 0.1796 ± 0.0010 siat ± 0.0084 sys , (32) 

where the systematic uncertainty was estimated to be 4.7% from the measurement on the 
branching ratio of K n2 in Section IIIIH8I It was noted that the trigger acceptance measured 
here was primarily due to geometry. 

Phase Space: The phase space acceptance was used to determine the acceptance of offline 
cuts on the momentum, range and energy. To measure the phase space acceptances (A PS ), 
events were first taken from those surviving from the trigger in the Monte Carlo simulation, 
and then the phase space cuts were applied. The ir + nuclear interaction was not included 
in this simulation. The acceptance was measured to be 

A PS = 0.3630 ± 0.0029 stat , (33) 

which included the loss due to both ir + absorption and decay-in-flight. 
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Correction for Nuclear Interaction: The nuclear interaction effect was investigated sepa- 
rately, in order to study the systematic uncertainty associated with it. K + — > 7i + uu events 
were generated with and without nuclear interaction, respectively. The ratio between the 
trigger acceptances multiplied by the ratio between the phase space acceptances gave the 
correction for nuclear interaction 

A nucl . = 0.4953 ± 0.0077 stat ± 0.0248^, (34) 

where the systematic uncertainty took into account the observed 0.15 cm difference on the 
range resolution (Table [TTJ) . This difference could affect the acceptance due to the Rdev 
cut and translated into a 5% uncertainty in the acceptance. It should be pointed that this 
definition took into account loss associated with nuclear interactions such extra energy in the 
detector associated with nuclear interactions that caused the PV counters to fire. This loss 
was not included in the i^-based acceptance Ar- 2 given in Table IXIl because muon-nuclear 
interactions are rare. 

6. Correction to T ■ 2 Trigger Inefficiency 

The T ■ 2 trigger required coincident hits both in the inner- most two RS layers and in the 
IC. It was noted that the T-2 simulation result did not include the acceptance loss due to the 
geometrical and counter inefficiencies of the T-Counters. The geometrical inefficiency was 
due to tracks passing through azimuthal gaps between adjacent T-Counters. The counter 
inefficiency occurred if the scintillation light induced by the charged track was not detected 
by the PMT's. This inefficiency was measured by using the and K^i events in the KB 
monitor data. UTC track extrapolation was required to give the expected T ■ 2 counter. 
The online T • 2 trigger condition of the reconstructed event data were checked to measure 
the T ■ 2 efficiency for and K n 2 events, separately. Since the energy losses in the T ■ 2 
counter from the and K n 2 events were different, simulations were done for these decay 
modes to obtain the average energy loss for and K n 2- Using an energy extrapolation 
gave the correction to the T • 2 inefficiency for the signal 

A T . 2 = 0.9358 ± O.OOll(stat) ± O.OUO(sys), (35) 

where the systematic uncertainty accounted for the fact that there was a 1.5% variation 
when changing the z requirement on the UTC track extrapolation to the T ■ 2 counter. 
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7. Normalization to the Branching Ratio 



Since a beam K + could decay after the Cerenkov counter with a daughter satisfying 
the B4 and target requirement in KB, or a beam K + could deposit energy in the B4 and 
target but exit the target without stopping, the total number of K +, s that satisfy the KB 
trigger requirement should be corrected for the K + stopping fraction (/ s ). This fraction was 
obtained by normalizing the total K + exposure to the K^ 2 branching ratio. The K^ 2 events 
were selected from the Kfi2 monitor trigger data with the same selection criteria as those 
used for signal, excluding all the cuts related to the ir + particle type, the BV and the BVL. 
The momentum, range and energy cuts for signal were replaced by a minimum 40 cm range 
requirement on the K^ 2 events. The K^ 2 acceptance measurement was also performed in 
the same way as that for the signal. This stopping fraction was computed by 

f = Nk^ 

Js N e K ff (K, 2 ) ■ Acc(K, 2 ) ■ B(K+ -> At +^) ' 

= 0.7740 ± 0.001 W, (36) 

where the 2 = 355, 119 and was the number of surviving K^ 2 events. N^^K^) = 
4.1475 x 10 6 and was the total exposure of K + 's {N^) with a prescaling factor for the 
Kfi2 monitor trigger during the data acquisition period. Acc^K^) = 17.4% and was the 
acceptance. 



8. Confirmation of the K n2 Branching Ratio 

Measurement of the K n2 branching ratio confirmed the validity of the evaluation of the 
acceptance of the K + — > ix + vu selection cuts. Stopped K n2 events were selected from the 
Kir2(l) monitor trigger data by imposing cuts similar to the K + — > tt + uu selection criteria 
except for those used for the PV cuts and for defining the kinematic signal region. Good 
K„ 2 events should also meet the requirements on the energy, momentum and range, which 
were defined to be within three standard deviations of the K v2 peaks. Fig. [321 shows the sta- 
bility of the measured K n2 branching ratio as a function of run number. The K n2 branching 
ratio was measured to be 

B (K + -> 7T+7T ) = — -j, ^ , 

N e K ff (K n2 )-Acc(K n2 )-f s 
= 0.219 ± 0.005 stat , (37) 
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FIG. 32: Measurement of the K+ — > 7r + 7r branching ratio as a function of run number in E949. 



where N^ n2 = 16, 405 and was the number of surviving 7^2 events. N e K ff {K n2 ) = 1.3233xl0 6 
and was the total exposure of K +, s {Nk) with a pre-scaling factor for the Kn2(l) monitor 
trigger during the data acquisition period. Acc{K^i) = 7.3% and was the acceptance. This 
branching ratio was in agreement with the world average [3j value of 0.209±0.001 within 
4.7%, indicating the validation of the acceptance measurement. This difference was treated 
as the systematic uncertainty and assigned to the acceptance A tr i g in Equation f l32l ). 



9. Summary of Acceptance and Sensitivity 

The acceptances of the K + — > -k + vv decay were split into several parts as given above. 
Table Kill summarizes all the contributions to the total acceptance in the standard region 
Atotai dard ■ To get the acceptance in the extended signal region, the estimated acceptance 
gains given in Section [HI G 71 were applied to the acceptances for A n ^^ e , Ak^ 2 and Ap$, 
yielding the final acceptance of 

Ace. = (2.22 ± 0.07 stat ± 0.15^,) x 1(T 3 . (38) 

This value is 10% higher than that in E787. It is noted that the acceptance of the standard 
E787 signal region for the E949 data was 84% of the acceptance for the E787 data due to 
losses incurred by the higher than expected instantaneous rates (Section III Ah . Based on 
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Contribution 


Acceptance 




0.3796 ± 0.0013 stat 




0.8785 ± 0.0029 stat 




0.6161 ± 0.0085 stat ± 0.0189 s . ys 


A 

-^n — >/i — >e 


0.3523 ± 0.0077 stat ± 0.0067 ays 


A-tvig 


0.1796 ± 0.0010 stat ± 0.0084^ 


Aps 


0.3630 ± 0.0029 stat 


A n ucl. 


0.4953 ± 0.0077 stat ± 0.0248 SJ/S 


At. 2 


0.9358 ± 0.0011 Btat ± 0.0U0 S y S 


fs 


0.7740 ± 0.001W 


Af o Z dard (xl0- 3 ) 


1.69 ± 0.05^ ±0.13 sya 



TABLE XII: A breakdown of the acceptance for the K + — > tt^vu selection criteria. 

the total exposure of K + 's (Nr), the single event sensitivity (SES) of the E949 2002 run 
was given by 

SES = (2.55 ± 0.08 stat ± 0.18^) x 10~ 10 . (39) 
I. Examining the Signal Region 

After the background analysis and the acceptance measurement were completed and 
satisfactory, all the selection criteria were then be applied to the data. At the stage of 
examining the signal region, no cut could be changed. 

One candidate event was observed inside the signal region. A close check also found 
that this candidate was located in the ir + — > fi + — > e + extended region as described in 
Section IIII G 71 Fig. [33] shows the range and kinetic energy of the events that passed all of 
the selection criteria, except for the phase space cuts on both the range and energy. This 
candidate together with the events observed in E787 are also shown in Fig. [33l As indicated 
in the Figure, the signal box definition in E949 was extended in comparison to that in E787. 
Fig. EH is an event display for this candidate. This event had a momentum of 227.3 MeV/c, 
a kinetic energy of 128.9 MeV and a range of 39.2 cm. Kinematically, this candidate event 
agreed with all the requirements as a signal, though it existed near the limit expected for 
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FIG. 33: Range vs kinetic energy of the events satisfying all of the cuts, except for the phase space 
cuts on both the range and energy. The plots are shown separately for E949 only (left) and E787 
plus E949 results. The rectangle represents the signal region defined in E787 (dashed lines) and 
E949 (solid lines). Events around E = 108 MeV were due to K n 2, which were not removed by the 
photon veto cuts. The light points in the right hand plot represent the expected distribution of 
K + — > it^vv events from simulation. 

signal. The measured quantities of the observed candidate used in the selection criteria 
were compared to the expected distributions for signal to evaluate the signal probability 
distributions for the candidate. The probabilities for the single beam K + requirements, the 
decay ir + kinematic requirements and the ir + — > fi + — > e + decay sequence cuts showed a 
fairly flat distribution, which was consistent with the expected signal distribution. There 
was no observed photon activity for this candidate. 



IV. RESULTS 



In this section we describe the method used to obtain the K + — > tx + vv branching ratio 
and the impact of the E949 and E787 K + — > 7i + uu candidates on the unitarity triangle. 



We also describe the implication of the results on the search for the 
K + — > tc + X° where X° is a stable, massless, non-interacting particle 



lypothetical decay 
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FIG. 34: Reconstruction of the candidate event (end view). The clusters of squares indicate both 
the K + track and the tt + track in the target. The hit IC sector is shown next to the ir + cluster in 
the target. The curve is the result of the UTC track fit. The circles along the track are the hits in 
the UTC. The radius of each circle gives the drift distance. The RS and RSSC hit layers are shown 
outside the UTC. Also displayed are the TD data in the ir + stopping counter, the reconstruction 
in the target and the CCD data in the K + stopping fiber. In the fit to the TD pulse shape, the ir + 
pulse (dashed) and the fi + pulse (dotted) are shown separately. No obvious 7r + pulse was observed 
in the CCD pulse shape for the K + stopping fiber. 

A. Background Functions 

We defined a number of cells in the extended signal region of differing signal/background, 
and calculated the expected signal/background (dA/dN) using individual A 1 ^ and N^fff 
functions for each rejection or background type. In total, we had seven types: TD rejection, 
PV rejection, K n 2 background, range tail background, /i + band background, single- 
beam background and double-beam background. The variation of each type for different cut 
position or cell could be expressed as the change of N^fff as a function of the corresponding 
acceptance A%*£, yielding seven functions in total used in this analysis. Two of them were the 
TD rejection vs. acceptance function (Fig. [26!) and the PV rejection vs. acceptance function 
(Fig. |29|) . Three of them were the relative background rate in the normalization branch vs. 
the relative acceptance functions for the kinematic background backgrounds (Fig. 1351) . The 
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FIG. 35: The expected relative kinematic background rate in the normalization branch vs. the 
relative acceptance for the K n 2 background (a), the range tail background (b) and the /i + 
band background (c). Both background rates and acceptance curves were normalized to one at the 
cut positions. 

rest were the relative background rate vs. the relative acceptance functions for the beam 
backgrounds (Fig. [361) . Both relative background rates and the relative acceptance curves 
were normalized to one at the cut positions. 



B. Likelihood Method 



The K + — > 7i + uu branching ratio was determined using likelihood analysis incorporating 
the predicted background functions in the signal region. The likelihood ratio X was defined 

as 

*=n ! — Ir^/-^- • < 40 > 

i=l 

where Sj and were the estimated signal and background in the i th cell, dj was the number 
of signal candidates in the i th cell and the product ran over all n cells 79J. In addition the 
likelihood estimator X f, s was defined as the value of X given the observed candidates. Cells 
were defined based on the predicted background functions described in previous section. 
The predicted background functions showed that there was additional background rejection 
capability within the signal region that could be exploited by sub-dividing the signal region. 
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FIG. 36: The expected relative beam background rate vs. the relative acceptance for the single- 
beam and CEX backgrounds (a) and the double-beam background (b). Both background rates 
and acceptance curves were normalized to one at the cut positions. 

The number of cells to be used for subsequent analysis, 3,781, was established prior to the 
examination of the signal region. 

The total background in the cell containing the signal candidate was estimated to be 
5.75 x 10 -5 dominated by a contamination of 4.92 x 10~ 5 events due to the K n2 range tail 
background. The ratio of the acceptance in this cell to the total acceptance in the standard 



region (A 



standard 
total 



0.001694) was estimated to be 1.21x10 4 . The expected number of 



signal events in this cell was 



N K -B- AfZt ard ■ M 



1.77 x 10 1Z x B x 0.001694 x 1.21 x 10' 
3.628 x 10 5 x B, 



(41) 



where B was the K + — > ti + vv branching ratio. 

C. Branching Ratio of K + — > tt + uv 
The central value of the branching ratio, denned as the value of B that maximized 

X bs 



80|, was 0.96 x 10~ 10 . Using only the E949 data, B(K + 



IX 1 vv\ 



(0.96j£$) x 10~ 10 



where the quoted 68% confidence level (CL) interval was determined from the behavior of 
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X as described in Ref. [79j and included only the statistical uncertainty. The estimated 
probability that the E949 candidate was due to background alone was 0.074. 

The results from E787 and E949 were combined to calculate the branching ratio for 
K + — > ir + uu. In the E787 K + — * ir + uu analysis, two K + — > 7i + uu candidate events were 
observed in the signal region j55|. The number of cells describing the E949 signal region 
were augmented by 488 cells that defined the signal region for the E787 analysis to produce 
a likelihood estimator X b s for the combined data. 

The confidence intervals for the combined E787 and E949 results took into account the 
estimated systematic uncertainties in the signal acceptance and the background rates. The 
systematic uncertainty of each background source was estimated to be about 15% based 
upon the results of the correlation studies. From the study in Section IHI H[ the systematic 
uncertainty on the acceptance was estimated to be about 8%. The systematic uncertainty 
of each background component and the acceptance were assumed to be uncorrelated and 
to follow a normal distribution with the magnitudes given above regarded as one standard 
deviation. With these assumptions, the K + — > 7i + uu branching ratio for the combined E787 
and E949 result was B(K + — > ix + vv) = (1.47^o|g) x 10~ 10 where the uncertainty denoted the 
68% CL interval. The corresponding 90% and 95% CL intervals were (0.27, 3.84) x 10~ 10 and 
(0.17,4.44) x 10~ 10 , respectively. The estimated probability that all the K + — > tt + uu can- 
didates observed in E787 and E949 were due to background was 0.001. The inclusion of 
the estimated systematic uncertainties had a negligible effect on the CL intervals due to the 
relatively poor statistical precision inherent in a sample of three candidate events. 



D. Search for B(K + ir+X°) 

The experimental signature of a K + — > tt + X° decay was identical to that of K + — > 7i + uu 
except that the kinematic signature afforded by the two-body decay (P^ = 227.1 MeV/c, 
E n = 127.0 MeV, R n = 38.6 cm) permitted the definition of a relatively high-acceptance, 
low-background signal region. The analysis method was identical to that used for the 
K + — > tt + uu analysis except that the signal region was defined to be within two standard 
deviations of the expected momentum, energy and range of the 7r + with the upper limits of 
the tightened to P < 229 MeV/c, E < 135 MeV, and R < 40 cm to suppress K + -> fi+X 
background. The expected background level was small (0.05 events), because the region was 
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far from the K n 2 peak. The acceptance studies for K + — > tt + X° decay paralleled those for 
the K + — > tt + uu. The single event sensitivity for the E949 K + — > tt + X° decay analysis was 
estimated to be (0.82 ± 0.02 stat ± 0.06 sys ) x lO" 10 . 

The candidate event observed in the signal region for K + — > n + vv E949 analysis was 
also in the K + — > tt + X° signal region. However, no candidates were observed in the 
K + — > 7T + X° signal region of E787 55J. The combined E787 and E949 sensitivity was 
0.196xl0~ 10 and, using the one observed candidate event without subtraction of the 
estimated background, the upper limit on the branching ratio was B(K + — > tt + X°) < 
0.73 x 10~ 10 at 90% CL using the Feldman-Cousins method 8jJ. This limit was larger 
than the previous 90% CL limit of 0.59 x 10~ 10 of E787 5fJ due to the E949 candidate 
event. 



E. Impact on the Unitarity Triangle 

As described in Section HJ the K + — > -k + vv branching ratio was directly related to the real 
and imaginary parts of \ t = V t * s V t d (Equation ([9])). In Fig. [37]the regions of the complex \ t 
plane allowed by the K + — > tt + vv branching ratio determined from the combined E787 and 
E949 results was compared to the regions allowed by other recent measurements with small 



theoretical uncertainties 82|. The region favored by other CKM-sensitive measurements is 
at the edge of the 68% CL region allowed by the K + — > tt + vv measurement. 

The other CKM-sensitive results 82J used to produce the confidence level intervals in 



Fig. [37] are dominated by measurements of B meson decays. The possible discrepancy 

B(K + —> 7i + vu) could 



83J , the clean theoretical 



between the Aj regions allowed by the S-decay measurements and 
be an indication of physics beyond the SM. As emphasized in Ref. 
interpretation of K — > iruu remains valid in most extensions of the SM in distinct contrast 
to the -B-decay measurements currently used to determine the CKM parameters. Thus a 
precise measurement of B(K + — > n + uu) would provide an unambiguous consistency test of 
the flavor sector of the SM. 
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FIG. 37: The allowed regions in the A( plane allowed by the combined E787 and E949 determination 
of the K + — > tt + uu branching ratio (gray), B + — > t + v (blue) and 5-mixing measurements (yellow). 
The regions outside the lighter (darker) shading have CL>0.95 (0.68). The area excluded by 
K + — > / k + uu at CL >0.95 is indicated by the gray line. The red-shaded region is allowed by the 
combination of these measurements and the small black region denotes the region allowed by all 



CKM-related measurements as evaluated by the CKMfitter Group 



V. CONCLUSION 



The rare decay K + — > -k + vv is a flavor-changing- neutral-current process and proceeds via 
1-loop diagrams mediated mainly by the top quark. Measuring B(K + — > tt + uu) is one of 
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the cleanest ways to extract \Vu\- 

In this paper we have reported results from the BNL experiment E949, an upgraded 
version of the BNL-E787 experiment, designed to improve the sensitivity for measurement 
of K + — > -k + vv decay. All the K + decays at rest were analyzed using a blind analysis 
technique in which the signal region was masked until the selection criteria were determined 
and the background levels were estimated. The development of the cuts and the estimation of 
the background levels were performed using a bifurcation method and, a likelihood analysis 
method was developed for interpreting the quality of candidate events. Enlargement of 
the signal region compared to E787 analysis increased the acceptance by 30% with a total 
background level in the signal region estimated to be 0.30 ± 0.03 events. 

An examination of the signal region yielded one event near the upper kinematic limit of 
the decay K + — > ir + ui?. Based on the candidate event, the branching ratio was determined 
to be B(K + tt+z/P) = (0.96+q;^) x 10~ 10 . E787 and E949 results were combined and 
the branching ratio was determined to be (1.47^Qgg) x 10~ 10 at the 68% CL level based 
on three events observed in the momentum region 211 < P < 229 MeV/c. The estimated 
probability that all the K + — > tt + uu candidates observed in E787 and E949 were due to 
background was 0.001. The measured branching ratio is in agreement with the SM prediction 
of (0.74 ± 0.20) x 10~ 10 within the uncertainty. 
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